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Abstract 

Context. Some circumstellar-interacting (CSI) supernovae (SNe) are produced by the explosions of massive stars that have lost mass 
shortly before the SN explosion. There is evidence that the precursors of some SNe Iln were luminous blue variable (LBV) stars. For 
a small number of CSI SNe, outbursts have been observed before the SN explosion. Eruptive events of massive stars are named as 
SN impostors (SN IMs) and whether they herald a forthcoming SN or not is still unclear. The large variety of observational properties 
of CSI SNe suggests the existence of other progenitors, such as red supergiant (RSG) stars with superwinds. Furthermore, the role of 
metallicity in the mass loss of CSI SN progenitors is still largely unexplored. 

Aims. Our goal is to gain insight on the nature of the progenitor stars of CSI SNe by studying their environments, in particular the 
metallicity at their locations. 

Methods. We obtain metallicity measurements at the location of 60 transients (including SNe Iln, SNe Ibn, and SN IMs), via 
emission-line diagnostic on optical spectra obtained at the Nordic Optical Telescope and through public archives. Metallicity val¬ 
ues from the literature complement our sample. We compare the metallicity distributions among the different CSI SN subtypes and 
to those of other core-collapse SN types. We also search for possible correlations between metallicity and CSI SN observational 
properties. 

Results. We find that SN IMs tend to occur in environments with lower metallicity than those of SNe Iln. Among SNe Iln, SN Iln- 
L(1998S-like) SNe show higher metallicities, similar to those of SNe IIL/R whereas long-lasting SNe Iln (I988Z-like) show lower 
metallicities, similar to those of SN IMs. The metallicity distribution of SNe Iln can be reproduced by combining the metallicity 
distributions of SN IMs (that may be produced by major outbursts of massive stars like LB Vs) and SNe IIP (produced by RSGs). 
The same applies to the distributions of the Normalized Cumulative Rank (NCR) values, which quantifies the SN association to Ft II 
regions. For SNe Iln, we find larger mass-loss rates and higher CSM velocities at higher metallicities. The luminosity increment in 
the optical bands during SN IM outbursts tend to be larger at higher metallicity, whereas the SN IM quiescent optical luminosities 
tend to be lower. 

Conclusions. The difference in metallicity between SNe Iln and SN IMs suggests that LB Vs are only one of the progenitor channels 
for SNe Iln, with 1988Z-like and I998S-like SNe possibly arising from LBVs and RSGs, respectively. Finally, even though line- 
driven winds likely do not primarily drive the late mass-loss of CSI SN progenitors, metallicity has some impact on the observational 
properties of these transients. 
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1. Introduction 


The study of the environment of supernovae (SNe) has become 
crucial to understand the link between different SN classes and 
their progenitor stars. Mass and metallicity are among the pro¬ 
genitor properties that can be investigated, and are fundamen¬ 
tal to understand stellar evolution and explosions (see |Anderson| 
|et al.|[20T5 for a review). The study of SN host galaxy metal¬ 
licity is now a popular line of investigation in the SN field. 
Metallicity can be obtained as a global measurement for a 


SN host galaxy (e.g. Prieto et al. 20081, following the known 
luminosity-metallicity or color-luminosity-metallicity relations 
(e.g., Tremonti et al. 2004 Sanders et al. |2013 1. It can also 
be estimated via strong line diagnostics when spectra of the 


* Based on observations performed at the Nordic Optical Telescope 
(Proposal numbers: P45-004, P49-016; PI: F. Taddia), La Palma, Spain. 


host galaxies are obtained. In particular, it has been shown 
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In this work we study the environments and in particular the 
metallicity of supernovae that interact with their circumstellar 
medium (CSM). A large variety of these transients has been ob¬ 
served in the past years, and in the following we briefly describe 
their physics, their observational properties and the possible pro¬ 
genitor scenarios, and how the study of their environments can 
help us to constrain the nature of their precursors. 
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Figure 1 CSI SN classification scheme. The SN explosion mech¬ 
anisms, the main elements in the spectra, the light curve proper¬ 
ties (only for SN Iln subtypes), the degree of association (NCR) 
to H II regions (only for SN IMs) and the possible progenitors 
are shown in correspondence with their CSI SN subtype. 

1.1. CSI SNe: subclassification and progenitor scenarios 

When the rapidly expanding SN ejecta reach the surround¬ 
ing CSM a forward shock forms which propagates through 
the CSM, while a reverse shock travels backward into the SN 
ejecta (see e.g., [Chevalier & Fransson||1994l l. The high energy 
(X-ray) photons produced in the shock region ionize the un¬ 
shocked CSM, giving rise to the characteristic narrow (full- 
width-at-half-maximumFWHM~100 km s“^) emission lines of 
Type Iln SNe (SNe Iln) (Schlegel |1990| l. Their emission lines 
are furthermore often characterized by additional components: 
a broad base (FWHM^IOOOO km s“^), produced by the ion- 


40000 km s-i) 
ized ejecta; and an intermediate (FWFIM~1000 km s“^) fea¬ 
ture, which originates in the cold- dense shell (CDS) between 
the forward and the reverse shocks. The kinetic energy of the 
ejecta is transformed into radiation, powering the luminosity of 
CSM-interacting (CSI) SNe. 

SNe Iln served as the first class of SN identified to be 
mainly powered by CSM interaction. Their CSM is H-rich, as 
revealed by their prominent Balmer emission lines. Within the 
SN Iln family, a wide range of observational properties is ob¬ 
served ( |Kiewe et al. 2012 Taddia et al. 2013a| l. Their optical 
light curves can decline sharply by several magnitudes (e.g., 
SN 1994W; |Sollerman et al.|199'8|l or settle onto a plateau lasting 
for years (e.g., SN 1988Z, Turatto et al.|1993[ SN 20 05ip, Fox 
et al.|2009||Smith et al.|2009b||Stritzinger et al.|2012| hereafter 
S12). 

A multiplicity of CSM geometries and densities, as well as 
a wide range of ejecta masses and kinetic energies (e.g., Moriya 
& Maeda||2014| l can explain this large variety of observational 
properties. This multiplicity may also suggest the existence of 
multiple progenitor channels for SNe Iln and other CSI SNe. 

The existence of subclasses within the large SN Iln group 
has recently been proposed ([Taddia et al.|2013a{ [Habergham et[ 


al. 


2014 hereafter H14). In the following paragraphs we con¬ 


sider SN Iln subtypes mainly based on the light curve properties, 
whereas other CSI SN types are defined based on their spectra. 


A scheme summarizing the different CSI SN types is shown in 

Fig- [3 

A group of long-lasting SNe Iln including SNe 1988Z, 
1995N, 2005ip, 2006jd and 2006qq was described in Taddia et 
[ar| ( [2013aj l. These are SNe whose light curves exhibit a slow de¬ 
cline rate of <0.7 mag (100 days)“^ over a period of time > 150 
days, or in some cases for several years. Their progenitors have 
been suggested to be Luminous Blue Variable stars (LBVs, see 
e.g. SN 2010jl, |Fransson et al.|2014| l. 

The LBV phase is traditionally interpreted as a transitional 
phase between an O-type star and a Wolf-Rayet (WR) star. 
According to standard stellar evolution models, a massive star 


should not end its life in the LBV phase. Recently, Smith & 
Tombleson (2015|l investigated the association of known LB Vs 


to O-type stars (which spend their short lives clustered in star¬ 
forming regions). They found that LB Vs are even more isolated 
than WR stars, and thus conclude that LB Vs are not single stars 
in a transitional phase between O-type and WR stars, but rather 
mass gainers in binary systems. In this scenario, the LBVs would 
gain mass from their massive companion stars which would then 
explode in a SN event that would kick the LBVs themselves far 
from their birth locations. When H14 and [Anderson et aLlP012| l 
compared the association of SNe Iln to bright Hll regions to 
that of the main CC SN types, it emerged that SNe Iln display 
a lower degree of association. As the degree of association to 
H II regions can be interpreted in terms of progenitor mass (see 
Sect. [7.3[ ), this result might suggest a low initial mass for the pro¬ 
genitors of SNe Iln, which apparently would be in conflict with a 
massive LBV progenitor scenario. The mechanism proposed by 
[Smith & Tombleson] ( [20151 ) reconciles the facts that LBVs are 
progenitors of at least a fraction of SNe Iln and at the same time 
they are found to have a weak positional association to bright 
H II region and O-type stars. 

There is indeed evidence that at least some SNe Iln arise 
from LBVs. Gal-Yam et al. ( 2007| l identified a likely LBV in pre¬ 
explosion images at the location of SN Iln 2005gl. SN 2009ip 
exhibited LBV-like outbursts before the last major event which 
might be explained with it s CC ([Prieto et aL 2013[ Mauerhan 
[et al.|[^13a[ [Smith et al.[[2014a| l. For several SNe Iln from 
the Palomar Transient Factory (PTF) there were signs of stel¬ 
lar outbursts before their terminal endpoint ( [Ofek et aL][2014| l. 
Furthermore, the CSM velocities observed in SNe Iln are of¬ 
ten consistent with those expected for LBV winds (~100- 
1000 km s“^), and the large mass-loss rates observed in SNe Iln 
are also compatible with those of large LBV eruptions, ranging 


between 10 ^ and 1 Mq yr ^ ( Kiewe et al.j2012 Taddia et al. 
2013a[ Moriya et al.|2014)l. The presence of bumps in the light 


curves of some SNe Iln can also be understood as the interaction 
with dense shells produced by episodic mass-loss events (e.g., 
SN 2006jd, S12). In summary, there is thus ample observational 
evidence from different lines of reasoning that connect LBVs 
and SNe Iln. 

However, also Red Supergiant stars (RSGs) with superwinds 
(e.g., VY Canis Majoris) have been suggested to be progeni¬ 
tors of long-lasting SNe Iln (see e.g., SN 1995N, Fransson et al. 
[2002[ Smith et al.|2009a| l, since they are characterized by strong 
mass loss forming the CSM needed to explain the prolonged in¬ 
teraction. 

There are other SNe Iln which fade faster than the afore¬ 
mentioned 1988Z-like SNe. These have early decline rates of 
~2.9 mag (100 days)“^, and given the linear decline in their 
light curves, these events can be labelled Type Iln-L SNe. 

Among the SNe Iln-L, it is possible to make further subclas- 
sifcations based on their spectra. Some SNe Iln-L show slowly 


2 











































































































F. Taddia et al.: Metallicity at the explosion sites of interacting transients 


evolving spectra (e.g., SNe 1999el, Carlo et al.|2()02l 1999eb, 


Pastorello et al. 20021, with narrow emission lines on top of 


broad bases lasting for months. Other SNe Iln-L display spec¬ 
tra with a faster evolution, with Type Iln-like spectra at early 
times and Type llL-like spectra (broad emission lines, some¬ 
times broad P-Cygni absorption) at later epochs (e.g., SN 1998S, 
|Fassiaetal.|2000l[20ml and SN 1996L, |B^ti et al.|2006| l. The 
different spectral evolution of these two subtypes can be traced 
back to the efficiency of the SN-CSM interaction and to different 
wind properties. Also for SNe Iln-L, both LBV (e.g., [Kiewe et| 
al.|2012|l and RSG (e.g., Frans son et al.|2005 1 progenitors have 


been considered in the literature. 

We note that both LB Vs and RSGs could be responsible for 
the production of CSI SNe when they are part of binary systems. 


We have already discussed the mechanism proposed by Smith & 
Tombleson ( 2015| l for LBVs. Mackey et al.| (2014l explain how 


interacting supernovae can come from RSGs in binary systems. 
In their scenario, the binary companion of a RSG can photoion- 
ize and confine up to 35% of the gas lost by the RSG during its 
life close to the RSG star itself, forming a dense CSM. When the 
RSG explodes, the SN will appear as a CSI SNe as the SN ejecta 
will interact with this dense shell. 

Mauerhan et al.|(2013b|l suggested the name “Iln-P” for ob¬ 


jects resembling SN 1994W, i.e., SNe with Type Iln spectra 
showing a ^100 days plateau in the light curve, similar to that of 
SNe IIP, followed by a sharp drop and a linear declining tail at 
low luminosity. The nature of these events is currently debated; 
some authors have proposed that 1994W-like SNe arise from the 
collision of shells ejected by the progenitor star and not from its 
terminal explosion (e.g., Humphreys et al.||2012[ Dessart et al. 


|2009| l. However, for SN 2009kn a bona-fide core-collapse (CC) 


origin was favoured by the observations of Kankare et a . (2012|l. 


sients could be something different from CSM interaction, such 
as radioactive decay of large amounts of ®®Ni or energy from 
a magnetar. We do not focus on the host galaxies of SLSNe II 
in this paper, with the exception of SN 2003ma, which shows 
a light-curve shape similar to that of SN 1988Z but brighter by 


^2.5 mag ([Rest et al.|201 Ijl. 

Besides SNe Iln, in the CSI SN grou 
called SNe Ibn, or 2006jc-like SNe (e.g.. 

3 we also find the so- 

Matheson et al.|2000 

Eoley et al. 

2007[ Pastorello et al.[[2007 

2008a[b[ Smith et al. 

2008a||Gor 

3ikov et al.|2014|l. These are transients showing He 


emission lines, arising from the SN interaction with a He-rich 
CSM. These SNe likely originate from the explosion of massive 
Wolf-Rayet (WR) stars. The WR progenitor SN 2006jc was ob¬ 
served to outburst two years before the SN explosion. 

Another class of objects arising from CSM interaction is 
the Type la-CSM or 2002ic-like subgroup. These events are in¬ 
terpreted as thermonuclear SNe interacting with H-rich CSM 


sum of narrow Balmer emission lines and SN la spectra diluted 
by a blue continuum ( Leloudas et al.|20T5 1 . 

Finally, a class of transients resembling SNe Iln is that of 
the SN impostors (SN IMs). These events are the result of out¬ 
bursts from massive stars. Most of the information regarding 
the observables of SN IMs has been collected by Smith et al. 
( |201 la[ see their table 9). In [Smith et al.| ( [TOl la| l SN IMs are 
interpreted as the result of LBV eruptions. Some LBVs may be 
observed in the S-Doradus variability phase (e.g. the 2009 op¬ 
tical transient in UGC 2773, [Smith et al.||2^ |2011a| l. LBVs 
in the S-Doradus phase are typically characterized by a vari¬ 
ability of 1-2 visual magnitudes, usually interpreted as the re¬ 
sult of a temperature variation at constant bolometric luminos¬ 
ity ( [Humphreys & Davidson) 1994} . Occasionally even giant 
eruptions of LBVs similar to that observed in Eta-Carinae dur¬ 
ing the 19th century may be observed in external galaxies (e.g. 
SN 2000ch, [Wagner et al.||20()4| [Pastorello et al.|[2010[ l. Their 
typical luminosities are lower than those of SNe Iln (M > —14), 
even though some bright events are now suspected to belong to 
this class (e.g. SN 2009ip, [Pastorello et al.[|2013[ [Fraser et al 


2013 Margutti et al.|2014[). Among SN impostors, objects like 


Others even suggested SN 1994W was the result of a fall-back 
SN creating a black hole ( [Sollerman|2002| l and it has also been 
proposed that SNe Iln-P might be electron-capture (EC) SNe 
( Nomoto|1984 l coming from the explosion of super asymptotic- 
giant-branch (AGB) stars (e.g., [Chugai et al.[[2004a[ [Mauerhan[ 
^t^I]20^. 

An additional member of the family of SNe Iln are the so- 
called superluminous SNe 11 (SLSNe II). These objects reach 
peak absolute magnitude of < —21 mag ( [Gal- Yarn 2012|l, and 
SN 2006gy serves as the prototypical example ( [Smith et al.[ 
2007^. We note that the mechanism powering these bright tran- 


SN 2008S have been proposed to be electron-capture SNe from 
super-AGB stars rather than LBV outbursts (e.g.,[Botticella et al. 
[2^ see Sect.[73]l. 

1.2. This work 

A possible approach to investigate if the LBVs, which produce 
(at least some) SN IMs, are the dominant progenitor channel for 
SNe Iln is to compare the properties of the environments for 
SNe Iln and SN IMs. In particular the metallicities of the sur¬ 
rounding material can be measured via strong emission-line di¬ 
agnostics. This is the focus of this paper. 

If the vast majority of SNe Iln originate from LBVs asso¬ 
ciated with SN IMs, then the metallicity distributions of these 
two groups should be similar. Else, if these distributions do not 
match, there could be room for other progenitor channels of 
SNe Iln. Metallicity measurements of SNe Iln and SN IMs based 
on the host-galaxy absolute magnitudes are provided by H14. 
In this paper we provide local metallicity estimates for a large 
sample of CSI transients, including SNe Iln, Ibn, la-CSM and 
SN IMs. Our measurements, carried out on data obtained at the 
Nordic Optical Telescope (NOT), are complemented with data 
available in the literature (e.g., Kelly & Kirshner[2012[ hereafter 


KK12, and H14) and in public archives. 

Metallicity measurements at the locations of CSI transients 
are important not only to compare the environments of SNe Iln 
and SN IMs, but also to clarify the role of metallicity in the 
prevalent mass losses of these events. It is well known that, given 
a certain stellar mass, a higher metallicity drives a larger mass 
loss in massive stars, due to stronger line-driven winds (e.g., 
[Kudritzki & Puls[[2000| l. However, the mass-loss rates due to 
this mechanism, observed in hot stars like Wolf-Rayet (WR), 
are on the order of 10“® Mq yr“^, which is substantially lower 
than what is required in SNe Iln. The massive CSM around 
SNe Iln must be produced by larger eruptions, whose underly¬ 
ing mechanism and metallicity dependence is largely unknown 
(e.g., gravity waves and super-Eddington winds have been in¬ 
voked, [Quataerr&^Shiodej20^|Smith]et^r|200^. Since the 
mass loss of SN Iln progenitors and SN IMs shapes their CSM 


(|Hamuy et al.||2003| Aldering et al.|2006 

IDilday et al.||2012| and thus the appearance of these transients, by looking for corre- 

[Taddia et al.[[2012[ [Silverman et al.[[2013[ 

Fox et al.[[2015|l, al- lations between metallicity and observational properties of these 

though a CC origin has also been proposed 

(Benetti et al.|2006 transients we aim to constrain to what extent metallicity is an 


important ingredient. In this paper we have collected the ob- 
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servables available in the literature for each SN with measured 
metallicity at its explosion location. We also complemented our 
analysis of the metallicity with that of another important prop¬ 
erty related to the SN environment, which is the association of 
SN locations to star-forming (SF) regions. In doing that, we used 
the results published by H14. 

The paper is organized as follows. In Sect. |^we introduce 
our SN sample and Sect, [^presents our observations and data 
reduction procedures. Section describes how we subtracted 
the underlying stellar population from each spectrum, how we 
measured the emission line fluxes and the spectral classification. 
Section concerns the method to obtain the local metallicity 
measurements, and Sect. presents the results on the metal- 
licities of CSI SNe, which are compared among the different 
CSI SN subtypes and with those of other CC SN classes. In 
Sect. 1^ we show the relations between metallicity and observ¬ 
ables of CSI SNe and finally the discussion and conclusions are 
given in Sect. and respectively. 


2. Sample of CSI transient host galaxies 

In Table[2we report the list of 60 transients included in our sam¬ 
ple. Thirty-five of them are SNe Iln, six are SNe Ibn, one is a 
SN la-CSM, 18 are SN IMs (if we count SN 2009ip also in the 
SN IMs, then we have 19 of these transients). 

With the NOT, long-slit spectra were obtained for the host 
galaxies of 13 SNe Iln, five SNe Ibn, one SN la-CSM, and 16 
SN IMs (the derived metallicity for SN 2007sv was published in 
Tartaglia et al.|201? i. The host galaxies observed at the NOT are 
marked with the letter “o” in the third column of Table [T] With 
the NOT we also obtained broad-band R and narrow-band Ha 
images for the SNe Iln (except for SN 1995G) and Ibn. 

The CSI transients observed at the NOT were chosen among 
those with published spectroscopic and photometric data. These 
SNe were thoroughly analyzed in the literature, and for most of 
them estimates of the physical properties of the CSM (e.g., wind 
velocity, mass-loss rate) were also available. This choice was 
made in order to select objects whose observational and phys¬ 
ical properties could be related to their host galaxy properties, 
such as the metallicity (see S ect. [Til. In par ticular, SN IMs were 
chosen among those listed by Smith et al. ( |2011a[ ). We also de¬ 
cided to observe nearby, i.e., resolved, host galaxies (z<0.026, 
see Table [g, to allow for the determination of the local metallic¬ 
ity. We discuss the possible biases introduced by this selection 
in Sect.m 

We complemented our observed sample with host galaxies 
whose metallicities (at the center of the galaxy or at the SN 
position) were already available in the literature (marked with 
the letter “1” in the third column of Table [^l or whose spec¬ 
tra were available in public archives (marked with the letter 
“a” in the third column of Table [T]). Eighteen SNe Iln (exclud¬ 
ing SNe 2005ip and 2006jd) have metallicities published in the 
literature (mainly from KK12 and HI4) and four have host- 
galaxy spectra obtained by the 6dF survey (Jones et al. 20091 and 
available via NEE|^ Eor two SNe Ibn their host-galaxy spectra 
are available in public archives. Metallicity measurements were 
available in the literature for nine SN IMs (excluding SN 2007sv 
that we have already published, and including SN 2009ip) and 
for six of them we also obtained observations at the NOT. Most 
of the metallicity estimates in the literature for SN IMs were ob¬ 
tained from P04. All the references are reported in Table (4th 
column). 


Nasa Extragalactic Database, http://ned.ipac.caltech.edu/ 


In summary, our entire sample includes a large fraction of 
the CSI SNe with published light curves and spectra. 

Considering the entire sample, on average SN IMs are lo¬ 
cated at lower redshift (< zjm > = 0.0020±0.0004) com¬ 
pared to SNe Iln (< zun > = 0.0276±0.0092) and Ibn (< 
zihn > = 0.0216±0.0078), since most of them were discovered 
only in nearby galaxies due to their lower luminosity. 


3. Observations and data reduction 

3.1. Nordic Optical Telescope observations and data 
reduction 


Observations at the NOT and data reductions were performed 
following the procedures outlined in |Taddia et al. (2013b i, where 
we studied the host galaxies of SN 1987A-like events arising 
from the explosions of blue supergiant stars (BSGs). In this sec¬ 
tion a brief summary of how our visual-wavelength spectroscopy 
was obtained and reduced is provided. 

Two observational campaigns (P45-004, P49-016) were car¬ 
ried out during 2012 and 2014. Six nights were spent observing 
the hosts of SNe Iln, Ibn and la-CSM in 2012 (four in April, one 
of which was lost to bad weather, and two in September). Four 
nights were spent observing the host galaxies of SN IMs (three in 
April 2014, one in September 2014). During April 2014, almost 
1.5 of the three nights were lost to bad weather. 

In both campaigns our main goal was to determine the metal¬ 
licity at the exact position of the CSI transients through strong 
emission-line diagnostics. In order to do that, we obtained long- 
exposure (> 1800 s), long-slit spectra of the Star forming (H ll) 
regions within the host galaxies, by simultaneously placing the 
slit at the SN position and through the galaxy center or through 
bright H II regions near the SN location. In order to intercept 
both the host galaxy center (or bright H II regions) and the SN 
location with the slit, we rotated the telescope field by the corre¬ 
sponding angle; then the slit was centered on a pre-determined 
reference star in the field and finally the telescope was offset to 
point to the SN location. The final pointing was checked with 
a through-slit image. In most cases the slit included a few H II 
regions, allowing for a determination of the host galaxy metal¬ 
licity gradient and hence of the metallicity at the distance of the 
SN from the host center (see Sect. [^. 

The instrumental setup that was chosen to acquire the host 
galaxy spectra at the NOT was the same as adopted for the study 
presented in [Taddia et al. p013b i, i.e. ALFOSC with grism #4 
(wide wavelength range ~3500—9000 A) and a 1.0"-wide slit 
(corresponding to the typical seeing on La Palma). The obtained 
spectral resolution is ~16-17 A. The exposure times adopted for 
each spectral observation are listed in Table 

The following procedure was adopted to carry out the spec¬ 
tral reductions. First, the 2D-spectra were bias subtracted and 
flat-held corrected in a standard way. When available, multiple 
exposures were then median-combined in order to remove any 
spikes produced by cosmic rays. We extracted and fitted with a 
low-order polynomial the trace of the brightest object in the 2D- 
spectrum (either the galaxy nucleus, or a bright star, or an H II 
region with a bright continuum). The precision of this trace was 
checked by plotting it over the 2D-spectrum. We then shifted 
the same trace in the spatial direction to match the position of 
each H II region visible in the 2D-spectrum, and then extracted 
a ID-spectrum for each H II region. The extraction regions were 
chosen by looking at the Ha flux profile, an example of which 
is presented in the top-left panel of Fig. where we also report 
the width of each spectral extraction. The extracted spectra were 
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Figure 2 (Top-right panel) Continuum-subtracted Ha image of UGC 4286, the host galaxy of SN 2010al. The 25th i?-band magni¬ 
tude elliptic contour is shown by a black solid line, along with the position of SN 2010al (marked by a red star) and the center of 
the galaxy (marked by a red circle). The slit position is shown, and a color code is used to present the N2 metallicity measurements 
at the position of each H II region that we inspected. (Top-left panel) Flux at the Ha wavelength along the slit, shown as a function 
of the distance from the SN location (marked by a dotted line, like the nucleus position). The N2 measurements are shown at the 
corresponding positions in the top sub-panel. (Bottom panel) Metallicity gradient of UGC 4286. The linear fit on our measurements 
is shown by a solid line. The interpolated metallicity at the SN distance is marked by a red square and its uncertainty corresponds to 
the fit error. The error bar (±0.2) for our N2 measurements is shown aside. The positions of SN and nucleus are marked by vertical 
dotted lines. The solar metallicity ( Asplund et al.|[2009 l and the LMC metallicity ( Russell & Dopita]|1990| l are indicated by two 
horizontal dotted lines. 


wavelength and flux calibrated using an arc-lamp spectrum and 
a spectrophotometric standard star observed the same night, re¬ 
spectively. Following [Stanishe^ ( |2007[ ), from each spectrum we 
removed the second order contamination, which characterizes 
the spectra obtained with grism #4. In this study we included 
all the spectra showing at least Ha, [N ll] A6584 and H/3 emis¬ 
sion lines. We identified the location of each corresponding H II 
region by inspecting the acquisition images. 

For SNe Iln, Ibn and la-CSM, we also had the chance to 
observe the host galaxies in a narrow Ha Alter and a broad R- 
band Alter, to build Ha maps as detailed in |Taddia et al.| ( 2013b| l. 
These maps further helped us to precisely determine the loca¬ 
tion of each H II region that was spectroscopically observed. An 
example of a continuum-subtracted Ha images is shown in the 


top-right panel of Fig. Here we indicate the SN position with a 
red star and the galaxy nucleus with a circle. An ellipse marks the 
25th i?-band magnitude contour. Each colored patch within the 
plotted slit aperture corresponds to the position of an extracted 
H II region spectrum. We also attempted to use these Ha maps 
to measure the degree of association (NCR, see Sect. |7.3[ ) to H II 
regions for each SN. However, since we had to cut the exposure 
times due to bad weather, we did not reach the desired depth 
in the Ha images. The time lost because of bad weather in the 
2014 campaign also prevented us from obtaining the same pho¬ 
tometric observations for the host galaxies of SN IMs. However, 
we used NCR data from H14 in the discussion about the SN Iln 
progenitor scenarios (see Sect. |8.1| l. Table summarize all the 
observations carried out during the two campaigns. 
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Rest wavelength [A] 

Figures (Top Panel) STARLIGHT stellar-population spectrum 
best-fit to the continuum of the spectrum obtained at the cen¬ 
ter of the host galaxy of SN 2006aa. The difference between 
the observed spectrum and the best fit gives the pure emission¬ 
line spectrum. (Bottom panel) The triple-Gaussian fit on the 
continuum-subtracted Ha and [N ll] lines of the bright H II re¬ 
gion at the center of the host galaxy of SN 2006aa. The observed 
fluxes are represented by black crosses, the best fit in red, the Ha 
component in blue and the [N ll] components in green. 


3.2. Archival data 

We used archival spectra to complement the dataset for SNe Iln 
and Ibn. All the four spectra retrieved via NED are low- 
resolution spectra obtained at the galaxy center by the 6dF 
Galaxy Survey ( |Jones et al.|2009| l. The spectrum of LSQ12btw 
was retrieved via the ESO archive and reduced. All these spectra 
included at least Ha, [N ll] A6584 and H/3 emission lines. 


4. Stellar population subtraction, line 

measurements and spectral classification 

In order to accurately measure the emission line fluxes necessary 
to determine the metallicity, the stellar population component of 
each spectrum was rn odelled with the STARLIGHT code ( |Cid| 
Fernandes et al. 2005 1 ^ and subtracted from each observed spec¬ 
trum (corrected for the Galactic extinction). This code linearly 
combines synthetic spectra from Bruzual & Chariot ( 2003| l and 
takes into account dust attenuation in order to fit the emission¬ 
line free regions of our spectra. The template spectra are stellar 
population spectra at hfteen different ages between 0.001 and 
13 Gyr and three metallicities (Z = 0.2, 1, 2.5 Zq). The same 
templates were used by Leloudas et al. (201 1[) in their study of 


http://astro.ufsc.br/starlight/ 


SNe Ib/c locations. As an example, the best-fit to the spectral 
continuum of SN 2006aa host galaxy nucleus is shown in Fig. 
(top-panel), as well as the “pure” emission line spectrum, as ob¬ 
tained from the observed spectrum subtracted of the best-ht con¬ 
tinuum. Once we obtained the pure emission line spectra, we 
proceeded to measure the fluxes of [N ll] A6584 and Ha, H/3 
and [O III] A5007. 

In order to measure the fluxes of the [N ll] A6584 and Ha 
emission lines, we followed a fitting procedure similar to that 
outlined in|Taddia et al.|P013bll. Given the low resolution of our 
ALFOSC spectra, we had to deblend [N ll] A6584, [N ll] A6548 
and Ha through a triple Gaussian fit. We fixed the width of each 
Gaussian to be the same, as determined by the spectral reso¬ 
lution. The known wavelength offsets between the centroids of 
the three lines was also fixed. The flux of [N ll] A6548 was fur¬ 
thermore fixed to be 1/3 of that in [N ll] A6584 (see |Osterbrock| 
& Ferland 20061. This assumption was needed to allow for a 


proper fit of this faint line, which could possibly contaminate 
the flux of Ha. Finally, the integrals of the two Gaussians fit¬ 
ted to [N II] A6584 and Ha provided us with the fluxes of these 
lines (see the bottom panel of Fig.[^for an example of the fitting 
procedure). 

For each spectrum we measured the H/3 and [O III] A5007 
line fluxes by fitting them with Gaussians. In several cases we 
could only place a limit on the flux of [O III] A5007. 

The ratios of line fluxes such as F([N ll] A6548) to F(Ha), 
and F([Olll] A5007) to F(H/3) can be used to determine the 
metallicity only if the dominant ionizing source for each re¬ 
gion are hot massive stars, since strong line diagnostics are 
based on this condition. To exclude other possible ionizing 
sources such as shock-excitation or AGN contamination, we 
classified our spectra based on the BPT diagram ([Baldwin"^ 


ffll981i, which shows logio(F([0 III] A5007)/F(H/3)) versus 


logio(F( N II] A6548)/F(Ha)), which is widely used in order 
to discriminate the excitation sources of emission line objects. 
logio(F([N II] A6584)/F(Ha)) is known as N2 and we will use 
it to determine the metallicity for our SNe (see Sect. |^. Gas 
ionized by different sources occupies different areas across the 


Sanchez et al. 

2014[ (Galbany et al. 20141; ( 

Belfiore et al. 

2015 

1). In Fig. 

4j we plot logio([0 III] A5007/H/3) versus N2 


for our spectra and checked if these points were located within 
the BPT diagram area given by |Kauffmann et al.| (|2003|l, i.e., 
logio([Olll] A5007/H^) < 0.61/(N2-0.05) -F 1.3, which de¬ 
fines the star-forming galaxies as opposed to galaxies with AGN 
contamination. The line ratios for each H II region are reported 
in Table We rejected the spectra falling in the other region 
of the BPT diagram. This occurred for the nuclear spectra of 
the host galaxies of SNe 1994W, 1995G, 2001ac, 2005dp and 
2005gl, and for another region of the hosts of SNe 200lac and 
1997bs. We also excluded the hosts of SNe 1996L (jBenetti et al. 


19991, 2005kj ( |Taddia et al.|2013^ and iPTFl liqb (Smith et al. 


20151 from our sample as we only had a nuclear spectrum with 
[Oiii] A5007/H/3) > 0.61/(N2-0.05) h- 1.3. 


5. Oxygen abundances at the SN explosion sites 

In this section we describe how the metallicity at the SN location 
was estimated from the observed and archival spectra after stellar 
population subtraction, line fitting and spectral classification. We 
also show how we included the metallicity values available in the 
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literature. For the observed spectra we followed the procedure 
illustrated in |Taddia et al.| ( |2013b| l. 


5. 1. Local metallicity measurements from observed and 
archival spectra 

Among the possible emission line diagnostics, we chose to use 
N2 ( Pettini & Pagel|2004l l for all our metallicity measurements. 
The oxygen abundance can be obtained from N2 using the 
following expression presented by Pettini & Pagel (|2004l: 


12+log(0/H) = 9.37+2.03 xN2+1.2xN22+0.32xN23. 

This expression is valid in the range —2.5 < N2 < —0.3 
which corresponds to 7.17 < 12+log(0/H) < 8 . 86 . This method 
has a systematic uncertainty of 0.18 which largely dominates 
over the error from the flux measurements. Following Thone et 
al. ( 2009| l and Taddia et al.| ( |2013'b l, we adopted 0.2 as the total 


error for each single measurement. We note that [Marino et al. 
( 2013| l have recently revised the N2 index using a large dataset of 
extragalactic H II regions with measured Tg-based abundances. 
However, we use the calibration by Pettini & Pagel ([2004]) to 


allow for a direct comparison with other SN types (see Sect. 6.2 


N2 has a larger systematic uncertainty than other methods 
such as 03N2 ( [Pettini & Pagel|2004| l. However, N2 presents sev¬ 
eral advantages that are briefly summarized here; 1 ) for nearby 
galaxies, the lines needed for this method fall in the part of 
the CCD with higher efficiency and thus are easier to detect 
than those needed for the other methods (e.g. 03N2 and R23, 
where lines in the blue or in the near ultraviolet are required). 
2) given the fact that Ha and [N ll] A6584 are very close in 
wavelength, this minimize the effects associated with differen¬ 
tial slit losses (we did not observe at the parallactic angle as we 
placed the slit along the SN location - host-galaxy center direc¬ 
tion) and uncertainties on the extinction. 3) It is well known that 
there are non-negligible offsets between different line diagnos- 


the metallicity for many other SNe in the literature (e.g.,|Thone 

et al. 

2009||Anderson et al.j 

20101|Leloudas et al.|2011|Sanders 

et al. 

2012 

Stoll etal.|2013| 

Kuncarayakti et al.[2013a|b Taddia 

et al. 

2013 

fl), this method al 

ows for a direct comparison of CSI 


transients to other SN types (see Sect. 

We note that there are possible drawbacks when using N2. 
For instance, the existence of a correlation between the metal 
licity derived from the N2 parameter and the N/O ratio ([Perez 


Montero & Contini|2009[l together with the presence of N/O ra¬ 


dial gradients across the disks of spiral galaxies (see e.g. P04; 
Molla et ^[2006 1 can affect the metallicity measurements ob¬ 


tained with this method. The most direct method to determine 
the oxygen abundances would be via the weak [O III] A4363 line 
flux (see e.g., Izotov et al.[[2006| l, but we could not detect this 
line in our spectra. 

For each Hll region with measured N2 metallicity and 
[Oiii] A5007/H/3) < 0.61/(N2-0.05) h- 1.3, and for each SN, 
we computed their de-projected and normalized distance from 
their host-galaxy nuclei, following the method illustrated by 
[Hakobyan etar] ( |2009l|2^ . 

In order to do that, we established the coordinates of the H II 
region/SN, and collected all the necessary information about its 
host galaxy; nucleus coordinates, major (2 R 25 ) and minor ( 2 b) 
axes, position angle (PA) and morphological t-type. All these 
data with the corresponding references are listed in Table [T] 
which also reports the result for the de-projected distance of each 
SN, rsAr/R 25 . We mainly used NED to collect the SN and galaxy 


coordinates as well as the galaxy dimensions and position an¬ 
gle. The Asiago Supernova Catalogue (ASC) was used for the 
t-type, SIMBACj^and HYPERLEDA^ was used when neither 
NED nor ASC included the aforementioned galaxy data. 

Having obtained the distance from the nucleus for each H II 
region (see the first column of Table [^, we plotted their metal- 
licities versus these distances, and performed a linear chi-square 
fit to the data. An example of such a fit is reported in the bot¬ 
tom panel of Eig. In most cases, the metallicity was found 
to decrease as the distance from the center increases (with the 
exception of nine hosts, where flat or positive gradients were de¬ 
rived). It is well known that there is such a negative metallicity 
gradient in galaxies (e.g., P04), likely because the star formation 
rate depends on the local density and therefore more metals have 
been produced in the denser, central parts of the galaxies. All the 
data and the fits are shown in Eigs. mi an d|7] 

The linear fits allowed us to interpolate or extrapolate the 
metallicity at the computed SN distance from the host center 
(see e.g. the red square in the bottom panel of Eig. |^. These val¬ 
ues were taken to be the local SN metallicity estimate. We note 
that these estimates were always found to match the metallicity 
of the H II region closest to the SN. In this way, we provided 
oxygen abundances also for those SNe that were not associated 
with a bright H II region (this is the case of several SN IMs, see 
Sect. [73]), and for those SN IMs or SNe Iln whose flux was still 
dominating the emission (e.g., SN IM 2000ch). All the local and 
central metallicities, and the metallicity gradients are reported in 
Table m 

When more than two H II regions were observed, we 
computed the uncertainty on the interpolated/extrapolated 
SN metallicity as the standard deviation of 1000 interpola¬ 
tion/extrapolations of 1000 Monte Carlo-simulated linear fits, 
based on the uncertainty on the slope and on the central metal¬ 
licity. These uncertainties are those rep orted in Table [4] and do 
not include the systematic N2 error (see Taddia et al.|2013b for 
a motivation of this choice). When only two measurements were 
obtained, we assumed 0.2 as the uncertainty for the metallicity 
at the SN position and at the host center. 

Eor most of the galaxies whose spectra were retrieved from 
online archives, only a single measurement was possible, often 
at the host center. In these cases, we adopted a standard metal¬ 
licity gradient (-0.47 R^ 5 ^) to extrapolate the metallicity at the 
SN position. This value corresponds to the average gradient of 
the large galaxy sample studied by P04. In these cases the er¬ 
ror on the SN metallicity included the extrapolation uncertainty, 
given by the ~ 0.1 RjTr;^ dispersion of the measured metallicity 
gradients (see Sect. [ 6 ^. 


5.2. Local metallicity values from the literature 

In the literature we found several published values for the metal¬ 
licities of CSI SN host galaxies. Some of them were local 
(e.g. a few from H14), most of them were measurements at 
the host galaxy center (e.g., all those obtained by KK12, ex¬ 
cept that of the host of SN 1994Y). KK12 presented 03N2 
metallicity estimates based on SDSS spectral measurement per¬ 
formed by an MPA-JHU collaboration and available online 
at http;//www.mpa-garching.mpg.de/SDSS/ Therefore we re¬ 
trieved the needed line fluxes for the N2 method for all the 
SN Iln hosts presented by KK12 and computed the metallicity 

^ https;//heasarc.gsfc.nasa.gov/W3Browse/all/asiagosn.html 
* http://simbad.u-strasbg.fr/simbad/ 

^ http: //atlas .obs -hp.fr/hyperleda/ 
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Figures (Top-left panel) Cumulative distribution functions (CDFs) of the central metallicities for SN Iln, Ibn and SN IM host 
galaxies. SN Ibn and SN IM hosts show slightly lower metallicities than SN Iln hosts. (Top-right panel) CDFs of the metallicity at 
the SN location for the same three classes. SN IMs and SNe Ibn are located in lower metallicity environments than those of SNe Iln. 
(Bottom panels) CDFs of the host metallicity gradients and de-projected SN distances from the host center for SNe Iln, Ibn, and 
SN IMs. 


measurements in the N2 scale. H14 and |Roming et al.] (2012i 
present five 03N2 values, so we had to convert them to the N2 
scale, using the relation presented by Kewley & Ellison| ( 2008| l. 
In the cases where the metallicity was measured at the host nu¬ 
cleus, we assumed the aforementioned standard gradient to ob¬ 
tain the metallicity at the SN de-projected distance. We note that 
we could have used the average gradient that we measured in 
some of our galaxies (see Sect. |^, instead of the average gradi¬ 


ent from P04. However, the sample of galaxies observed by P04 
is larger and representative of many morphological types, so we 
prefer to adopt their value. This assumption could potentially af¬ 
fect the results concerning SNe Iln (for 13 of the 35 SNe Iln we 
adopt the P04 average gradient), whereas it is not important for 
most of the SN IMs or SN Ibn, where we do not need to assume 
any average gradient. However, we will show in the following 
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Figure 9 CDFs of the metallicity at the SN location for CSI transients and for other SN classes from the literature. The legend entries 
have been ordered by mean metallicity. 


sections that the SN Iln results are solid independently of the 
average gradient assumption. 

6. Metallicity resuits 

In the following sections we describe the metallicity results 
for our sample of CSI transients, including comparisons among 
CSI SN subtypes and to other CC SNe. 

6.1. Metallicity comparison among CSI transients 

In Table we report the average values for the metallicity at the 
host center, at the SN position, and for the metallicity gradient. 
The mean gradients are computed excluding the objects where 
the gradient was assumed to be -0.47 R^ 5 ^. 

SNe Iln show slightly higher average oxygen abundance 
at their locations (8.47±0.04) than do SN IMs (8.33±0.06). 
The latter values are consistent with the metallicity of the 
Large Magellanic Cloud (LMC has 8.37, [Russell & Do pita 
1990 1 , which is sub-solar (the Sun has 8.69, [Asplund et al. 
2009 I. To understand if this difference is statistically signifi¬ 
cant, computing the mean values is not sufficient. Therefore, 
we need to compare their cumulative distributions, which are 
shown in the top-right panel of Fig. through the Kolmogorov- 
Smirnov (K-S) test. We found that the probability that the 
metallicity distributions of SNe Iln and SN IMs are differ¬ 
ent by chance is only 1%. If we exclude the CSI SN hosts 
where we assumed the average P04 metallicity gradient to de¬ 
termine the local oxygen abundance, this difference still holds 
(p-value = 1.5%), with < 12 -flog( 0 /H)> 77 „ = 8.48±0.04 and 
< 12-Hog( 0/H) >jM = 8 31±0 .07. Also if we use the N2 cali¬ 
bration by [Marino et al.j ([2013[), the K-S test gives a p-value of 


1%. To account for the uncertainties of each metallicity measure¬ 
ment (see Table|^, we performed the K-S test between 10® pairs 
of Monte Carlo-simulated metallicity distributions of SNe Iln 
and SN IMs. The 68 % of the resulting p-values were found to 
be <0.17, suggesting that the difference is real. We discuss this 
result and its implications for the SN Iln progenitor scenario in 
Sect.[^ SNe Ibn show an average metallicity (8.44±0.11) simi¬ 
lar to that of SNe Iln and a distribution that appears in between 
those of SNe Iln and SN IMs, although here our sample is lim¬ 
ited to only 6 SNe Ibn and there is no statistically-significant 
difference with the other two distributions (p-values>56%). 

The average oxygen abundances at the host center for 
SNe Iln and SN IMs are 8.63±0.03 and 8.59±0.06, respectively. 
These distributions (top-left panel of Fig. turned out to be sim¬ 
ilar, with p-value=47%. SN Ibn hosts exhibit a slightly lower 
mean central metallicity than SNe Iln (8.55±0.11, but with p- 
value=54%). Metallicities at the center are expected to be higher 
than those at the SN explosion site due to the negative metallicity 
gradient typically observed in the galaxies (e.g., P04). 

The mean gradient for the host galaxies of SNe Iln, 
Ibn and SN IMs is -0.10±0.09 R 25 \ -0.33±0.09 R^g^ 
and —0.31 ±0.06 R^s^, respectively. These absolute values are 
slightly lower than the average gradient (—0.47 R^g^) of the 
galaxies studied by P04 and adopted for the SNe in our sample 
with a single metallicity estimate at the galaxy center. The gra¬ 
dient distribution of SN Iln hosts shows higher values than those 
of SN IMs and SNe Ibn (bottom-left panel of Fig. [^, however 
the K-S tests result in p-values > 12%. 

SNe Iln, Ibn and SN IMs show similar average distances 
from the host center, with 0.64±0.15 R 25 , 0.43±0.10 R 25 and 
0.54±0.08 R 25 , respectively. Their distributions are almost iden¬ 
tical (p-values>49%, see the bottom-right panel of Fig. [^. 
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Figure 10 (Left-hand panel) SN Iln light curves from the literature. Three main subclasses are identified (plus the SLSNe) based 
on the light curve shape. Long-lasting SNe Iln (1988Z-hke; blue) show slow decline (~0.7 mag (100 days)“^) and sustained 
luminosity for >5 months. SN Iln-P (1994W-hke; red), show a month plateau followed by a sharp drop, similar to what 
is observed in SNe IIP. SNe Iln-L (1998S or 1999el-hke, depending on the spectral evolution; magenta) show a faster decline 
(^2.9 mag (100 days)“^). References for the light curve data and the extinction of SNe with measured metallicity can be found in 
Table[5] For SNe 1994aj, 1999el, 2005cl, 2006gy, 2006tf and 2009kn the references are Benetti et al. ( |1998[ l, Di Carlo et al. ( |2002| l, 
Kiewe k al.|P012l l, Smith et al. (20071, [Smith et al. ( |2008b| l, and Kankare et al. (2012 i, respectively. (Right-handpanel) CDFs of 


the metallicity at the SN location for our CSI transients, including three SN Iln subclasses (1988Z-hke, 1994W-hke, 1998S-hke). 
1998S-hke SNe are located at higher metallicity compared to the other subtypes. 


We note that if we had assumed our average metallicity 
gradient instead of that from P04 for the hosts with a single 
measurement at the center, then the metallicity distribution of 
SNe Iln (which is the only one that would be affected by a differ¬ 
ent choice of the average metallicity gradient) would have been 
at slightly higher metallicity, making the difference with SN IMs 
even more significant. 


7. Host galaxy properties and CSI transient 
observables 

In the following sections we investigate any potential correla¬ 
tions between the observables of CSI transients and the proper¬ 
ties of their environment, in particular the metallicity. This exer¬ 
cise aims to establish to what extent metallicity plays an impor¬ 
tant role in shaping the appearance of these events. 


6.2. Metallicity comparison to other SN types 


7.1. Metallicity and SN properties 


Figure shows the cumulative distributions of the metallicity 
at the SN location for our CSI transients, along with those of 
other SN types previously published in Taddia et al. (2013b, 
see hgure 15 and references in the text) and [Leloudas et al.| 
( 2014| l. All the metallicities included in the plot are local and 
obtained with the N2 method in order to enable a reliable com¬ 
parison among the different SN classes. It is evident that the 
SN IMs have relatively low metallicities along with SN 1987A- 
like events, SNe Ilb and SLSN 11. When we compare the metal- 
licity distribution of SN IMs to that of SNe IIP, their difference 
has high statistical significance (p-value = 0.09%; even higher 
if the comparison is to SNe Ibc). SNe Iln show higher metal- 
licities than SN IMs, but still lower than those of SNe IIP (p- 
value=29%) and SNe Ibc (p-value=2.5%, highly signihcant). 
All the SN types show higher metallicities compared to SLSNe I 
and SLSNe R. 


7.1.1 . Metallicity and SN Iln subtypes 

In the introduction we described three SN Iln subgroups, 
mainly based on the light curve shapes, i.e., 1988Z-hke (or 
long-lasting SNe Iln), 1994W-hke (or SNe Iln-P), and 1998S- 
like (or SNe Iln-L with fast spectral evolution). A collection 
of SN Iln light curves from the literature, including most 
of the light curves of SNe Iln with measured metallicity, is 
shown in Fig. (left-hand panel), where we display the 
different subtypes with different colors. In our sample we 
selected the events belonging to the different SN Iln subclasses 
(see the second columns of Tables m to look for possible 
trends in the metallicity. In Fig. (right-hand panel) we plot 
their metallicity CDF. It was possible to subclassify only 21 
of the 35 SNe Iln in our sample, i.e. those with sufficient 
photometric coverage in the literature. We found that our 
SN Iln-L (we only have 1998S-hke SNe) are typically located at 
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solar metallicity (<12+log(0/H)>77„_i = 8.69±0.03). 

Long-lasting SNe lln are found at lower metallicity 
(<12-Flog(0/H)>77n-88Z = 8.38±0.06), similarly to SN IMs. 
Among them, the object showing the largest oxygen abundance 
is SN 2005ip. The difference between SN 1998S-like events and 
long-lasting SNe lln is statistically significant (p-value = 0.004) 
and would be even larger if we had adopted our measured 
average metallicity gradient instead of that from P04, since we 
assumed the P04 gradient for the hosts of three (out of five) 
SNe Iln-L. It would also be statistically significant if we exclude 
the hosts where we assumed the P04 gradient from both groups. 
1994W-like SNe exhibit a metallicity distribution similar to that 
of 1988Z-like SNe «12-Flog(0/H)>77n-p = 8.39±0.12). In 
this group we include only SNe 1994W, 2006bo and 201 Iht, 
which do show a plateau followed by a sharp decay. 

7.1.2. Metallicity and SN Iln/Ibn magnitude at peak 

It has been noticed in the literature that some bright SNe Iln 
exploded in low-metallicity environments (Stoll et al. ]20 TT] i. In 
this paper, we have the opportunity to test if such a trend can be 
confirmed for a larger sample of SNe Iln. For each CSI SN we 
collected their peak apparent magnitudes in five optical bands 
{Ulu, B, V, i?/r/unfiltered, I/i) when available in the litera¬ 
ture. We also coll ected the visual Galactic extin ction for each SN 
(Ay {MW), from Schlafly & Finkbeiner 2011 1 and, when avail¬ 


able in the literature, an estimate of the host extinction { Av{h)). 

Peak magnitudes and extinctions are summarized in Table ^ 
with references provided. We computed the distance from the 
redshift (see Table [^, assuming Hq = 73.8 km s“^ Mpc“^ 
( |Riess et al.||201l] l. Using the distance and the extinction es¬ 
timate, from the apparent peak magnitude of each filter we 
computed the absolute peak magnitude {Muiu,B,v,Rlr junj.,i 
which are reported in Tablej^and plotted against the local metal¬ 
licity estimates in Fig. [m When the absolute peak magnitude is 
only an upper limit, it is marked with a triangle. This figure does 
not show any clear correlation or trend between metallicity and 
luminosity at peak. When we group our events in different metal¬ 
licity bins (i.e., log(0/H)+12 < 8.2, 8.2 <log(0/H)H-12 > 8.4, 
8.4 < log(0/H)4-12 > 8.6, log(0/H)+12> 8.6) and compare 
their absolute magnitude distributions via K—S tests, no statisti¬ 
cally significant difference is found. We thus see no evidence that 
normal SNe Iln exhibit brighter absolute magnitudes at lower 
metallicities. However, we note that SLSNe II (M^/tj < -21), 
which might also be powered by CSI, are typically found at 
lower N2 metallicities (8.18 < log(0/H)-)-12 > 8.62) than those 
of our normal SN Iln sample ( [Leloudas et al.|20T4) l. 

7.1.3. Metallicity and dust emission 


|Fox et ak] ( |201 l| l observed a large sample of SNe Iln at late 
epochs in the mid-infrared (MIR) with Spitzer. Among their 
SNe Iln, 17 events are included in our sample and six of them 
show dust emission (including SN 2008J that was later re-typed 
as a la-CSM; |Taddia et al. 2012[). We compare the metallicity 


faint SN. However, a too high radiation field may also destroy the 
circumstellar dust. Indeed we see a statistically-significant dif¬ 
ference in r/i?/Mn/.-band peak absolute magnitudes between 
the SNe with and without dust emission (p-value = 5%), with 
the dust-rich SNe being brighter by 1.3 mag on average, (see 
Fig.[^ right-hand panel). Here we include the objects with only 
a limit on the peak magnitude in the computation of the CDFs. 

7.1.4. Metallicity and SN Iln/Ibn CSM properties 

For several of our SNe Iln and Ibn, we collected data about their 
CSM/progenitor wind velocities (v^,, as measured from the blue- 
velocity-at-zero-intensity BVZI, FWHM or P-Cygni absorption 
minima of their narrow lines), and their mass-loss rates {M, es¬ 
timated from spectral analysis, e.g. Kiewe et al.|2012| or via light 
curve modelling, e.g. |Moriya et aL|2014[ ). We list these quanti¬ 
ties in Table|^ with their references. In Fig.[T^(left-hand panels) 
we plot M and Vu, against the local metallicity measurements. It 
is important to stress that the mass-loss rates are sometimes ob¬ 
tained by measuring the shock velocity from the broad compo¬ 
nent of the emission lines. This might be misleading in some 
cases, as the broadening is sometimes produced by Thomson 
scattering in a dense CSM ( Fransson et al.|20T4l l. 

If we consider only mass-loss rates >10“^ Mq yr“^, a trend 
of higher M with higher metallicities is visible in the top-left 
panel of Fig. This behavior is also what is expected in hot 
stars, where line-driven winds drive the mass loss, but this is 
typically for much lower mass-loss rates. The metallicity depen¬ 
dence of the mass-loss rates due to line-driven winds is a power- 
law (PL) with a = 0.69 (e.g., Vink 201 l| l. We show this PL with 
a black segment in Fig. and we can see that it is roughly 
consistent with the observed in the data. When we compare 
the mass-loss rate CDFs of metal-poor and metal-rich SNe Iln, 
the difference between the two populations is found to be sta¬ 
tistically significant (p-value = 0.030-0.034, depending on the 
metallicity cut, see the top-right panel of Fig. [T^. We found five 
SNe with relatively low mass-loss rate (M < 10“^ Mq yr“^) 
that seem to deviate from the aforementioned trend. These are 
SNe 1998S and 2008fq, which belong to the SN Iln-L subclass, 
and SNe 2005ip and 1995N, two long-lasting (i.e., 1988Z-like) 
events. We note that all the SNe IIL whose metallicity measure¬ 
ments are available in the literature on the N2 scale ( [Anderson et| 
ari|2010 l show relatively high metallicity (log(0/H)-Fl2~8.6), 
similar to those of SNe Iln-L. Furthermore, the SNe IIL class 
is characterized by relatively weak CSM interaction and rela¬ 
tively low mass-loss rates up to 10“^ Mq yr“^ (e.g., SN 1979C 
and SN 1980K, Lundqvist & Fransson|1988] l, again resembling 
SNe Iln-L, which typically show lower mass-loss rates than 
SN 1988Z-like events (see the top-left panel of Fig. [T2| . In 
this sense, SNe Iln-L might be considered as close relatives to 
SNe IIL (SN 1998S was sometimes defined as a SN IIL, Chugaij 
& Utrobin 2008| l, having only slightly stronger CSM interac¬ 
tion, as suggested by the persistent emission lines in their early 
phases. 


distribution of those with dust emission against those without 
dust emission, and do not find any statistically significant differ¬ 
ence (see Fig. 13 left-hand panel), despite the fact that one could 


expect the formation of more dust at higher metallicity. However, 
we notice that the detection of dust emission is likely more af¬ 
fected by other factors, such as the luminosity of the SN itself. 
Since the optical SN radiation is reprocessed and emitted in the 
MIR, a luminous SN can reveal the presence of dust better than a 


In the bottom-left panel of Fig. 12 there are evidence that 
SNe Iln at higher metallicity tend to have larger wind velocities 
than those located in low metallicity environments. However, 
when we compare the wind velocity CDFs of metal-poor and 
metal-rich SNe Iln (see the bottom-right panel of Fig. [I^ , this 
difference is only marginally significant (p-value = 0.14-0.27). 
SNe Ibn show higher wind velocities than SNe Iln, regardless 
the metallicity. We note that, in line-driven winds, the metallic¬ 
ity dependence of the wind velocity is weak and can be described 


11 
















































F. Taddia et al.: Metallicity at the explosion sites of interacting transients 


by a PL with a = 0.12 (|Kudritzki|2002| l. We plot this PL with 
a black segment in Fig. 12 and we can see that is roughly con¬ 
sistent with the observed trend. SNe Iln-L (98S-like) typically 
show higher wind velocities than long-lasting SN Iln ( 88 Z-like). 

We discuss these results and their implications on the mass- 
loss mechanism in Sect.f 


7.2. Metallicity and SN IM properties 


The comprehensive paper by [Smith et al.| ( |201 la| l lists the ab¬ 
solute magnitudes of the outburst peak(s), the absolute magni¬ 
tude of each SN IM progenitor before the outburst, the expan¬ 
sion velocity during the outburst (Wexpi from spectral analysis), 
the characteristic decay time (ti. 5 , i.e. the time the SN IM spent 
to fade by 1.5 mag from peak), the Ha equivalent width (EW), 
the presence of multiple peaks and sharp dips in the light curves. 
[Smith et al.| ( [201 la| l noted that there are no clear correlations 
among these observables. 

On the other hand, we have found some indications of possi¬ 
ble trends with respect to the metallicity at the SN IM location. 
In Fig. [^(central panel), the SN IM progenitor absolute magni¬ 
tude before the outburst is plotted versus the metallicity. Squares 
correspond to detections, triangles to upper limits. The differ¬ 
ent colors correspond to the different pass-bands (see caption). 
It is evident how (in the optical) most of the fainter SN IM pro¬ 
genitors are located at higher metallicity, and the brighter ones 
at lower metallicity, with the exception of SN 1961V (and of 
those events where we only have poor limits, > —12 mag). Note 
that [Smith et al. ( 2011a[) suggest a CC origin for this event. On 
the other hand, [Van Dyk & Matheson[ ( [2012[ ) suggests that the 
1961V event was the outburst of a LBV progenitor star as it is 
still detectable in recent HST images. 

In the bottom panel of Fig. the difference in magnitude 
(AM) between the progenitor of the outburst and the peak of the 
eruptive event is shown versus the metallicity. The higher the 
metallicity at the location of the transient, the larger seems to be 
the difference between the peak luminosity during the outburst 
and the luminosity of the progenitor. 

As both the SN IM progenitor magnitudes and AM show 
a possible trend with the metallicity, it follows that these two 
observables also show a trend, with lower AM values for the 
more luminous SN IM progenitor. 

The large variations in optical magnitudes (AM > 2) ob¬ 
served in these SN IMs are probably due to actual variations in 
the bolometric luminosities, rather than being an effect of varia¬ 
tions in temperature at constant bolometric luminosity as in most 
of the Galactic LB Vs, which are characterized by smaller AM 
(see e.g., Clark et al.|2009| and reference therein). However, we 
also note that AM might be affected, or even produced, by the 
formation of dust in the CSM surrounding the quiescent SN IM 
progenitor ( |Kochanek|20 1 4|l. 

The top panel of Fig. |l4| reveals that there is no clear trend 
between peak outburst luminosity and metallicity. The same is 
true for expansion velocity, Ha EW and ti 5. Furthermore, we 
found no difference in the metallicity distributions of SN IMs 
with and without multi-peaks in the light curves, or between 
SN IMs with and without sharp dips in the light curves. 


7.3. Association to SF regions and CSI transient 
properties 

The degree of association of a SN to star forming (SF) region; 
in a galaxy can be interpreted in terms of SN progenitor zert 



Figure 15 CDFs of the metallicity at the SN location for SN IM 
subtypes, mainly based on NCR values. 


age main sequence mass (James & Anderson|2006[ Anderson & 


[James|2d08| [Anderson et al.|2012| l, i.e., the stronger the associa 

tion, the larger the progenitor mass. A way to quantify the degree 
of association is to measure the NCR index, i.e. the normalized 
cumulative rank pixel value. This number is obtained by ranking 
the pixels in the Ha continuum-subtracted host-galaxy image, 
based on their counts. Then, the CDF of the pixel counts is built 
and the NCR value is the CDF value of the pixel corresponding 
to the SN position. 

H14 presented a comprehensive study of the association to 
SFRs for CSI transients, where they introduced a subclassihca- 
tion of the SN IMs; SN 2008S-like and p Carinae-like. The hrst 
ones (SNe 2008S, 1999bw, 2001ac, 2002bu, 2010dn, 2008-OT) 
do not show any association to bright Hit regions (NCR=0) 
and are often enshrouded in dusty environments. The second 
category (SNe 1954J, 1961V, 1997bs, 2002kg and VI) are as¬ 
sociated with SF regions (NCR>0). Also because of the weak 
association to H II regions, the hrst class of SN IMs has been 


suggested to be related to 8-10 Mq super-AGBs (Botticella et 
ar|[2009[ Pumo et al. 2009 [ Wanajo et al. 2009 Prieto et al. 
20i0[ [Kochanek[[201 1 1 rather than to more massive LBV pro 


genitors (however, see Smith & Tombleson 2015| l, which are 
more likely to produce 77 Carinae-like IMs. When comparing the 
metallicity distributions of these two populations, we do not hnd 
any statistically signihcant difference (see Fig.[T5]l. However the 
numbers are small and the uncertainties are large. The same is 
true for SNe Iln and Ibn. If we plot the SN IM properties from 


Smith et al. (2011ai (see previous section) against the NCR val¬ 


ues from H14, we do not hnd any clear trends or correlations. 
Furthermore, the NCR does not exhibit any clear trends with 
peak absolute magnitudes, mass-loss rate or CSM/wind velocity 
in SNe Iln. In the last column of Table[4|we list all the available 
NCR values from H14. 
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We note that |Crowther| ( |2013| l concluded that the observed 
association between CC SNe and Hit regions provides only 
weak constraints upon the progenitor masses. 


8. Discussion 


8.1. Implications for the SNI In progenitor scenario 

As described in the introduction, there are several indications 
that LB Vs, which likely produce SN IMs, are also progenitors 
of at least a fraction of SNe Iln. On the other hand, the large 
variety of SNe Iln may suggest the existence of multiple precur¬ 
sor channels. The latter is also suggested by our result on the 
different metallicity distributions of SNe Iln and SN IMs. Since 
we know that at least some SNe Iln are produced by the termi¬ 
nal explosion of LBVs (see e.g., Pastorello et al.|2007 Gal-Yam 
et al. 2007 Ofek et al. [2014 1 , and assuming that most of the 
SN IMs are produced by LBVs, then the metallicity distribution 
of SNe Iln can be considered to be a combination of the metal¬ 
licity distribution of SN IMs with that of one or more alterna¬ 
tive progenitor populations. These other progenitors must have 
higher metallicity than that of SN IMs, which would explain why 
SNe Iln are found at higher metallicity compared to SN IMs. 

RSG stars with superwinds like VY Canis Majoris have been 
invoked as possible progenitors of some SNe Iln (e.g., Fransson 
|et al.|[2002l [Smith et al.||2009a] ). If all the SNe Iln which are 
not produced by LBVs are actually produced by RSGs, and if 
we assume that the metallicity distribution of SNe IIP (which 
do have RSG progenitors, Smartt||2009 1 and that of RSGs pro¬ 
ducing SNe Iln are the same, then we can combine the metal¬ 
licity distributions of SN IMs and SNe IIP to reproduce the ob¬ 
served metallicity distribution of SNe Iln. This exercise will al¬ 
low us to roughly estimate what fractions of LBVs and RSGs 
that are precursors of SNe Iln, in a simple two-progenitor pop¬ 
ulations scenario. This is done by maximizing the p-value of 
the K-S test between the SN Iln metallicity distribution and that 
of the combined distributions of SN IMs/LBVs-fSNe IIP/RSGs. 
To produce combined metallicity distributions of SN IMs/LBVs 
and SNe IIP/RSGs, we linearly interpolate their observed CDFs 
(which are marked by a green dashed line and a red solid line in 
the top-left panel of Fig.[T6| varying the number of objects in the 
two samples, while always keeping a total of 35 events, which 
is the number of SNe Iln that we have in the metallicity CDF 
(solid black line in the top-left panel of Fig. [T6| . This is done 
to produce a meaningful comparison to the sample of SNe Iln. 
The two extreme cases are obtained when the combined distri¬ 
bution contains only SN IMs/LBVs or only SNe IIP/RSGs. In 
between, we obtain 33 different distributions with mixed frac¬ 
tions of SNe IIP/RSGs and SN IMs/LBVs. We then performed 
the K-S test between all these distributions and that for SNe Iln, 
and plot the resulting p-values as a function of the SN IMs/LBV 
fraction in the top-right panel of Fig. 
p-value reaches its maximum between 
SN IMs/LBVs as SN Iln progenitors, with the rest being formed 
by SN IIP/RSGs. This range would change to 20%-35% if we 
exclude the hosts where we adopted the gradient from P04. The 
CDF formed by 40% of SN IMs and 60% of SNe IIP/RSGs is 
shown with a dashed line in the top-left panel of Fig. [T^ and it 
is a good match to the distribution of SNe Iln. We notice that, 
among SNe Iln, those resembling SN 1998S show a metallicity 
distribution similar to those of SNe IIP and IIL, suggesting that 
this SN Iln subclass might preferentially arise from RSGs rather 
than from LBVs. On the contrary, 1988Z-like SNe show a metal- 
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It is evident that the 
28% and -48% of 


licity distribution consistent with that of the SN IMs, suggesting 
a LBV origin. 

By using data from H14 and Anderson et al.| ( 2012| l, a simi¬ 
lar exercise was repeated with the NCR distributions for SNe Iln, 
SNIMs and SNe IIP. Wefound(see the bottom panels of Fig. [T6| 
that the SN Iln NCR distribution can be reproduced with high 
statistical signihcance by a simulated population of SN IMs 
(from LBVs, -40%) and SNe IIP (from RSGs, -60%), iden¬ 
tical to what was found when comparing the metallicity distri¬ 
butions. It is interesting to note that jOfek et al. ( |2014| l estimate 
that —50% of SNe Iln show outbursts only a few years before 
collapse, which is similar to our fraction of SN IMs/LBVs. This 
might imply that most LBV progenitors of SNe Iln display out¬ 


bursts just before collapse. Indeed it has been found by Moriya 
et al. P014|l that SN Iln progenitors suffer enhanced mass loss 


as they get closer to collapse. 

In reality, the scenario is likely more complicated, espe¬ 
cially if we consider that some SN IMs, i.e. the SN 2008S-like 
events, are likely produced by super-AGB stars, and not by LBV 
eruptions. However, the metallicity distribution of SN IMs like 
SN 2008S is very similar to that of SN IMs like rj Carinae (which 
are likely produced by LBV eruptions). In this perspective, we 
can consider the estimated 40% of LBV progenitors as an up¬ 
per limit. It is more realistic to think that of this 40% some 
SN Iln have for instance super-AGB stars as precursors. SNe Iln- 
P (1994W-like events), which indeed show a metallicity distri¬ 
bution similar to SN IMs (see Fig. left-hand panel), might be 
the result of EC-SNe from super-AGB stars. 


If we look at SN rates, SN II (II P and IIL) and SNe Iln 
form —55% and —9% of the CC SNe ( [Smith et al.|2011bl l, re¬ 
spectively. In our simplihed two-progenitor population scenario, 
—60% of the SNe Iln comes from RSGs with large mass-loss 
rates (similar to the massive hypergiant VY Canis Majoris). If 
we assume a Salpeter initial mass function and single stars as 
progenitors, to reproduce the ratio between SNe Iln from RSGs 
and SNe II (—10%) would require that RSG stars with initial 
mass between 8 and 20 Mq produce SNe IIP/L and massive 
RSG stars with initial mass between 20 and 25 Mq produce 
SNe Iln. These ranges are reasonable, but they do not account for 
binary systems, which might be important to explain the progen¬ 
itors of SNe Iln ( [Smith & Tombleson|2015| l and other CC SNe. 

There are potential biases in the samples of SNe Iln and 
SN IMs, which might affect the comparison between their metal¬ 
licity distributions. Most of the SNe in our samples were discov¬ 
ered by galaxy-targeted surveys, with the exceptions of a few 
SDSS (Sloan Digital Sky Survey) II SN survey and PTF SNe. 
Targeted surveys often pitch larger galaxies with more metal rich 
distributions that the non-targeted surveys ( [Sanders et al.|20f^ 
see e.g.,), so both these distributions are likely biased in that di¬ 
rection. 


The SNe Iln are also located at larger redshifts than the 
SN IMs. The latter are more difficult to discover at large dis¬ 
tances given that they are typically —4 magnitudes fainter than 
SNe Iln. For the redshifts in this investigation, there is however 
no metallicity evolution expected (see e.g. Savaglio 2006[ l. There 
is also no clear correlation between the SN Iln luminosity and 
their host metallicity, with the brighter SNe Iln (more likely to 
be discovered at larger distances) showing similar host metallic¬ 
ity to the fainter ones. 


As most SNe Iln are selected based on the availability of lit¬ 
erature data, it could be that the SN Iln sample is biased towards 
peculiar objects, which are more often discussed in the litera¬ 
ture as compared to the “normal” ones. It is, however, unclear 
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Figure 16 (Top-left panel) Metallicity CDFs of SNe Iln, SN IMs, SNe IIP and of a simulated SN progenitor population constructed 
of SN IM/LBVs and SN IIP/RSGs, which fit the observed SN Iln metallicity CDF. (Top-right panel) p-value of the K-S test between 
the metallicity distributions of SNe Iln and the simulated progenitor population as a function of the fraction of LBV s in the simulated 
progenitor population. (Bottom-left panel) NCR CDFs of SNe Iln, SN IMs, SNe IIP and of a simulated SN progenitor population 
made of SN IM/LBVs and SN IIP/RSGs, which fit the observed SN Iln NCR CDF. Data from H14 and [Anderson et aT] ( |20I21 l. 
(Bottom-right panel) p-value of the K-S test between the NCR distributions of SNe Iln and the simulated progenitor population as 
a function of the fraction of SN IM/LBVs in the simulated progenitor population. 


if and how this selection would affect the observed metallicity 
distribution. 

An alternative interpretation of the difference between the 
metallicity distribution of SN IMs and SNe Iln might be the fol¬ 
lowing. If we instead assume that all the SNe Iln are actually 
produced by LBV stars, then the SN IMs are preferentially pro¬ 
duced by the LBV at low metallicity, whereas at relatively high 
metallicity LBV would not produce giant outbursts but only ex¬ 
hibit e.g., S-Doradus variability. However, we do know that large 
eruptions like p Car and other LB Vs are found in both the MW 
and other nearby metal-rich galaxies (e.g., M 31, M 81, M 101), 
as well as in nearby metal-poor galaxies such as NGC 2366 and 


the SMC (see |Weis|2006| for a list of LB Vs and their metallic- 
ities). We therefore do not favour this explanation for the lower 
metallicity of SN IMs as compared to that of SNe Iln. 

8.2. Implications for the mass loss mechanism 

The large mass-loss rates derived for SN Iln progenitors and ob¬ 
served in SN IMs are not explicable solely in terms of steady, 
metal line-driven winds ( |Smith|2014b| l. This mechanism cannot 
reproduce M > 10“^ Mq yr“ , which we typically deduce for 
these transients. Eruptive mechanisms must play an important 
role. 
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The relation between metallicity and outbursts has not yet 
been clarified, but here we present a number of related obser¬ 
vational results. We need also to consider that, via steady line- 
driven winds, the metallicity might anyway have an effect on the 
appearance of CSl transients, as they act for the entire life of 
their progenitors. 

Despite the large spread, in SNe Iln it is possible to see that 
higher mass-loss rates and faster CSM velocities correlate with 
higher metallicity at the SN location (Fig.[T^. Note that most of 
the mass-loss rate estimates used assume steady winds as well as 
shock velocities derived from the width of the broad component 
of Ha, even if the progenitor winds may in fact not be steady 
(Dwarkadas 20111, and the broad component of the emission 
lines can be due to Thomson scattering in a dense CSM (|Chugai| 
jet al.|2002| l. 

The bottom panel of Fig. shows that in SN IMs, a higher 
metallicity tend to prompt relatively brighter outbursts in the op¬ 
tical. AM might be produced by an actual variation in the bolo- 
metric luminosity of the SN IM or be the effect of a temperature 
variation at almost constant bolometric luminosity as is observed 
for “S Doradus” LB Vs (e.g., [Humphreys & Davidson|1994||Vink 


2011 1 . These LBVs present temperatures of several lO'^ K in the 


quiescent state, and then become F-type stars (with temperatures 
of 8-9 X 10^ K) during the outburst. In this case, at larger metal¬ 
licity we would have found the SN IMs whose photospheric tem¬ 
peratures suffered the stronger cooling during the outburst, i. e., 
the SN IMs whose progenitors had the higher temperatures be¬ 
fore the outburst. More bolometric information on our SN IMs 
are necessary to draw further conclusions. 

The fact that we found relatively brighter outbursts at larger 
metallicity might fit with the mass-loss bi-stability jump sce¬ 
nario discussed by |Vink| ( |201 Ij ) for LBVs. When a LBV becomes 
cooler than 2.5 x lO'* K, the Fe recombines from Fe IV to Fe III, 
enhancing the opacity and thus the mass loss. LBVs are likely to 
cross the bi-stability temperature threshold several times during 
their lives, inducing variable mass loss. Obviously, large metal- 
licities would favor this mass-loss enhancement and hence the 
luminosity in the eruptive phase. 

Figure (central panel) shows that progenitors of SN IM 
outbursts at lower metallicity tend to have higher optical lumi¬ 
nosity. If we assume similar temperatures for the SN IM progen¬ 
itors, this would mean that those with higher bolometric lumi¬ 
nosity tend to be more metal poor. However, if we instead as¬ 
sume that they have similar bolometric luminosities, the brighter 
the progenitor in the optical, the lower would its temperature be. 
Therefore a possible implication of this result is that the cooler 
progenitors of SN IMs tend to have lower metallicities. Again, 
more data are needed to establish the bolometric and color prop¬ 
erties of SN IMs. 

The fact that at higher metallicity SN IMs show larger AM 
but their progenitor show lower absolute magnitudes could also 
be explained if we take into account the presence of dust in 
SN IMs. A larger metallicity would favour larger mass loss and 
thus stronger CSM interaction (and larger AM), but at the same 
time more dust would surround the star and thus the progeni¬ 
tor would suffer of larger extinction. However, there are SN IMs 
with large AM both with (SN 2008S, SN 2010dn, OT2008) and 
without (SNe 1997bs and 2003gm) dust emission ([Thompson et 
|al.|2009D . 


9. Conclusions 

We summarize our main findings as follows; 


• SNe Iln are located at higher metallicities than SN IMs, and 
this difference is statistically significant. 

• The locations of SNe Iln-L (1998S-like) exhibit the highest 
metallicities among SNe Iln. Their metallicity distribution is 
similar to those of SNe IIL and IIP (produced by RSG pro¬ 
genitors). On the other hand, long-lasting SNe Iln (1988Z- 
like) are typically metal-poorer and exhibit a metallicity dis¬ 
tribution similar to that of SN IMs (that may be produced by 
LBV progenitors). 

• The metallicity distribution of SNe Iln can be interpreted 
as a combination of the metallicity distributions of SN IMs 
and SNe IIP, which might be similar to the metallicity dis¬ 
tributions of LBVs and RSG stars, respectively. The same is 
true for the NCR distribution of SNe Iln. These results on 
the metallicity and NCR distributions are consistent with a 
scenario where SNe Iln have both LBV (~40%) and RSG 
(^60%) progenitors. If we also consider the possibility that 
some SN IMs come from super-AGB stars, then the esti¬ 
mated ~40% fraction of LBV progenitors might be consid¬ 
ered as an upper limit. 

• We do not find a significant difference in the metallicity dis¬ 
tributions of CSI SNe with and without dust emission in the 
MIR, but rather a difference in the distributions of the peak 
absolute magnitude, with the dust-rich SNe being brighter. 

• Above M ~ 10“^ M© yr“^, SNe Iln located at higher 
metallicities tend to show larger mass-loss rates. There are 
also indications suggesting larger CSM/wind velocities for 
SNe Iln at higher metallicity. 

• There are hints that at higher(lower) metallicity the differ¬ 
ence in optical magnitudes between the SN IM progenitor 
and its outburst’s peak tend to be larger(smaller) and SN IM 
progenitors tend to have lower(higher) optical luminosity. 
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Table 1. Sample and basic data of 63 (38 SNe Iln, 6 SNe Ibn, 1 SN la-CSM, 18 SN IMs) CSI SN host galaxies. 
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2007SV IM I/O UGC5979 0.003723 10:52:40.05 +67:59:14.2 10:52:41.16 +67:59:18.8 1.6983 1.0190 9.9'^ 80 NED HYP 0.18 

2008S IM I/O NGC6946 0.000133 20:34:45.35 +60:05:57.8 20:34:52.32 +60:09:14.1 11.4817 9.7594 5.9 53 NED NED(IR) 0.63 

2008- OT IM 1 NGC300 0.000480 00:54:34.16 -37:38:28.6 00:54:53.48 -37:41:03.8 21.8783 15.5336 7.0 111 NED NED 0.43 

2009- OT IM o UGC2773 0.000771 03:32:07.24 +47:47:39.6 03:32:07.32 447:47:38.6 1.7000 1.3090 10.0 105 NED(B) NED(IR) 0.03 

2010da IM 1 NGC300 0.000480 00:55:04.86 -37:41:43.7 00:54:53.48 -37:41:03.8 21.8783 15.5336 7.0 111 NED NED 0.21 
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F. Taddia et al.: Metallicity at the explosion sites of interacting transients 


Table 2. Log of the photometric and spectroscopic observations at the Nordic Optical Telescope with ALFOSC. 


SN Type Galaxy Photometry Spectroscopy 

Filter Integration (sec) Integration (sec) 


1987F 

Iln-88Z 

NGC 4615 

R 

3x40 

1x1800 




Ha(#50) 

3x180 


1988Z 

lln-88Z 

MGC+03-28-022 

R 

3x30 

1x1800 




Ha(#51) 

3x150 


1994W 

lln-P(94W) 

NGC 4041 

R 

3x40 

1x1800 




Ha(#49) 

3x200 


1995G 

lln-88Z 

NGC 1643 



2x1800 

1995N 

lln-88Z 

MGC-02-38-016 

R 

3x40 

1x1800 




Ha(#49) 

3x200 


1998S 

lln-L(98S) 

NGC 3877 

R 

3x40 

1x1800 




Ha(#21) 

3x150 


2005db 

Iln 

NGC 214 

R 

3x60 

3x1800 




Ha(#50) 

3x240 


2005cp 

Iln-88Z 

UGC 12886 

R 

3x90 

3x1800 




Ha(#51) 

3x300 


2005gl 

Iln-L(98S) 

NGC 266 

R 

3x90 

3x1800 




Ha(#50) 

3x300 


2005ip 

lln-88Z 

NGC 2906 

R 

3x60 

1X 3600 




Ha(#49) 

3x200 


2006aa 

Iln 

NGC 3947 

R 

3x40 

2x1800 




Ha(#51) 

3x180 


2006bo 

Iln-P(94W) 

UGC 11578 

R 

3x60 

3x1800 




Ha(#50) 

3x240 


2006jd 

Iln-88Z 

UGC 4179 

R 

3x40 

1200+1800+1452 




Ha(#50) 

3x200 


2008J 

la-CSM 

MGC-02-07-033 

R 

3x60 

3x1800 




Ha(#50) 

3x240 


1999cq 

Ibn 

UGC 11268 

R 

3x40 

1x1800 




Ha(#51) 

3x200 


2002ao 

Ibn 

UGC 09299 

R 

3x40 

2x1800 




Ha(#49) 

3x180 


20051a 

Ibn/Iln 

KUG 1249+278 

R 

3x40 

2x1800 




Ha(#50) 

3x200 


2006jc 

Ibn 

UGC 4904 

R 

3x40 

2x1800 




Ha(#49) 

3x200 


2010al 

Ibn 

UGC 4286 

R 

3x40 

2x1800 




Ha(#50) 

3x200 


1954J 

IM 

NGC 2403 



1300 

1961V 

IM 

NGC 1058 



3x1800 

VI 

IM 

NGC 2366 



1800 

1997bs 

IM 

NGC 3627 



1800+900 

1999bw 

IM 

NGC 3198 



1800 

2000ch 

IM 

NGC 3432 



1800+2400 

200 lac 

IM 

NGC 3504 



1800+900 

2002bu 

IM 

NGC 4242 



1800 

2002kg 

IM 

NGC 2403 



1300 

2003gm 

IM 

NGC 5334 



1800+1363 

OT2005 

IM 

NGC 4656 



1800 

2006fp 

IM 

UGC 12182 



3x1800 

2007SV 

IM 

UGC 5979 



1800 

2008S 

IM 

NGC 6946 



3x1800 

2009-OT 

IM 

UGC 2773 



4x1800 

2010dn 

IM 

NGC 3184 



1800+900 


Note. — Spectroscopic observations were performed with Slit 1.0^^+Grism #4. We could not image the host 
galaxies of the SN IMs because during the observational campaign of April 2014 several nights were lost due 
to bad weather, and we decided to spend the remaining nights doing spectroscopy. 
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Table 4. 


Metallicity estimates and Normalised Cumulative Rank (NCR) pixel index for the sample of CSI transients. 


SN 

Type 

Data 

source 

Reference 

12+log(0/H) 
(N2 at the 
host center) 

12-blog(0/H) 
(N2 at rsisr) 

12+log(0/H) 

gradient 

NCR 

1987F 

Iln-88Z 

0 

This work 

8.67(0.07) 

8.47(0.18) 

-0.27(0.21) 

0.352 

1988Z 

IIn-88Z 

o/a 

This work + SDSS 

8.47(0.01) 

8.34(0.05) 

-0.23(0.08) 


1994W 

lIn-P(94W) 

0 

This work 

8.60(0.05) 

8.61(0.08) 

+0.05(0.21) 

0.795 

1994Y 

lIn-L(98S) 

1 

KK12 


8.74(0.20) 

-0.47* 

0.000 

1994ak 

Iln 

1 

H14^ 


8.74(0.20) 


0.000 

1995G 

Iln-88Z 

0 

This work 

8.39(0.11) 

8.46(0.14) 

+0.14(0.15) 


1995N 

IIn-88Z 

o 

This work 


8.31(0.18)* 


SN 

1997ab 

Iln-88Z 

1 

KK12 

8.35(0.18) 

8.20(0.20) 

-0.47* 


1998S 

lIn-L(98S) 

o/a 

This work + SDSS 

8.62(0.07) 

8.67(0.13) 

0.06(0.14) 


1999gb 

lln 

1 

H14^ 


8.66(0.20) 


0.676 

2002A 

lln 

1 

H14^ 

8.82(0.20) 

8.66(0.20) 

-0.47* 

0.401 

2002ea 

Iln 

1 

KK12 

8.78(0.18) 

8.65(0.20) 

-0.47* 


2003ke 

Iln 

1 

KK12 

8.81(0.18) 

8.58(0.21) 

-0.47* 


20031O 

Iln 

1 

H14^ 


8.61(0.20) 


0.000 

2003ma 

IIn-88Z 

1 

iStoll etal.((2011( 


8.45(0.18) 



2004F 

Iln 

1 

"KKrz —"—' 

8.71(0.20) 

8.60(0.20) 

-0.16 


2005bx 

IIn-L(98S) 

a 

NED 

8.78(0.18) 

8.67(0.20) 

-0.47* 


2005cp 

IIn-88Z 

0 

This work 

8.60(0.03) 

8.53(0.05) 

-0.51(0.24) 


2005db 

Iln 

o 

This work 

8.35(0.05) 

8.51(0.06) 

+0.45(0.08) 

0.398 

2005gl 

IIn-L(98S) 

0 

This work 

8.53(0.20) 

8.59(0.20) 

+0.17 

0.000 

2005ip 

IIn-88Z 

l/o 

S12 


8.73(0.18)“ 


SN 

2006aa 

Iln 

o 

This work 

8.71(0.08) 

8.75(0.14) 

+0.12(0.30) 


2006am 

Iln 

1 

KK12 

8.50(0.18) 

8.28(0.21) 

-0.47* 

0.000 

2006bo 

IIn-P(94W) 

0 

This work 

>8.69 

<8.37 

< -0.65 


2006db 

Iln 

1 

KK12 

8.33(0.18) 

8.02(0.21) 

-0.47* 


2006jd 

IIn-88Z 

l/o 

S12 


8.28(0.18)“ 


SN 

2006qq 

IIn-88Z 

a 

NED 

8.78(0.18) 

8.08(0.25) 

-0.47* 


2007rt 

IIn-88Z 

1 

KK12 

8.82(0.18) 

8.74(0.20) 

-0.47* 


2007va 

Iln 

1 

iKozlowski et al. 1(2010) 


<8.08 



2008B 

Iln 

1 

'kK12 -'*-' 

8.67(0.18) 

8.46(0.21) 

-0.47* 


2008fq 

IIn-L(98S) 

a 

NED 

8.80(0.18) 

8.76(0.20) 

-0.47* 


2009ip 

IM-^IIn 

1 


Margutti et al.|j2014j 

8.80(0.02) 

8.25(0.10) 

-0.36(0.07) 


2010jl 

IIn-88Z 

1 


Stoll etal.li2Un) ' 


8.25(0.18) 



201Iht 

IIn-P(94W) 

1 


Roming et al.|(2012|'= 


8.20(0.18) 



2013L 

IIn-88Z 

a 

T?ED 

8.68(0.18) 

8.17(0.23) 

-0.47* 


2008J 

la-CSM 

o 

This work 

8.76(0.07) 

8.73(0.07) 

-0.60(0.25) 

0.807 

1999cq 

Ibn 

0 

This work 

8.88(0.05) 

8.81(0.05) 

-0.48(0.10) 


2002ao 

Ibn 

o/a 

This work 

8.38(0.02) 

8.21(0.06) 

-0.40(0.12) 


20051a 

Ibn/IIn 

o 

This work + SDSS 

8.56(0.09) 

8.18(0.25) 

-0.46(0.28) 


2006jc 

Ibn 

0 

This work 

8.26(0.20) 

8.25(0.20) 

-0.01 


2010al 

Ibn 

o 

This work 

8.65(0.02) 

8.57(0.03) 

-0.32(0.07) 


LSQ12btw 

Ibn 

a 

ESO 


8.62(0.18) 



1954J 

IM 

l/o 

P04 

8.52(0.11) 

8.42(0.11) 

-0.35 

0.187 

1961V 

IM 

l/o 

P04 

8.71(0.04) 

8.43(0.04) 

-0.32 

0.363 

VI 

IM 

o 

This work 


7.84(0.18)t 


0.000 

1997bs 

IM 

o 

This work 


8.66(0.21) 

-0.47* 

0.023 

1999bw 

IM 

l/o 

P04 

8.69(0.07) 

8.35(0.07) 

-0.64 

0.000 

2000ch 

IM 

o/a 

This work + SDSS 

8.34(0.07) 

8.24(0.12) 

-0.09(0.08) 

SN 

200lac 

IM 

0 

This work 

8.77(0.20) 

8.79(0.20) 

0.02 

0.000 

2002bu 

IM 

o 

This work 

8.49(0.20) 

8.29(0.20) 

-0.24 

0.000 

2002kg 

IM 

l/o 

P04 

8.52(0.11) 

8.45(0.11) 

-0.35 

0.055 

2003gm 

IM 

0 

This work 


8.39(0.21) 

-0.47* 

0.000 

2005-OT 

IM 

o 

This work 


<7.76** 



2006fp 

IM 

0 

This work 


8.18(0.21) 

-0.47* 

0.965 

2007SV 

IM 

l/o 


Tartaglia et al.||2015| 


8.05(0.05)tt 



2008S 

IM 

l/o 


PW -- 

8.70(0.06) 

8.44(0.06) 

-0.41 

0.000 

2008-OT 

IM 

1 

P04 

8.49(0.13) 

8.32(0.13) 

-0.40 


2009-OT 

IM 

o 

This work 

8.19(0.04) 

8.20(0.04) 

0.11(0.10) 


2010da 

IM 

1 

P04 

8.49(0.13) 

8.40(0.13) 

-0.40 


2010dn 

IM 

l/o 

P04 

8.97(0.06) 

8.75(0.06) 

-0.63 

0.000 

<IIn> 


0 + 1 +a 


8.63±0.03 

8.47±0.04 

-0.10±0.09 


<Ibn> 


o+l+a 


8.55±0.11 

8.44±0.11 

-0.33±0.09 


<IM> 


o+l+a 


8.59±0.06 

8.33±0.06 

-0.31±0.06 



* The metallicity gradient is assumed to be -0.47 R 25 ^, from the average of the metallicity gradients presented by P04. 

** Metallicity from the N2 limit obtained from the spectrum of the closest Hll region, as the gradient of the galaxy is difficult to 
obtain given that the host of SN 2005-OT is strongly titled and irregular. 

^Metallicity obtained from the spectrum of the closest HII region, as the gradient of this large galaxy is difficult to obtain given that 
the host of VI is larger than the ALFOSC field of view. 

“Metallicity from the spectrum of the closest Hll region as measured in SI2, as we could not determine the emission line fluxes 
from our NOT spectra, where the SN is still detected. 

^ ^Metallicity from the spectrum of the closest H11 region, see |Tartaglia et al.|2015| 

^Metallicity from the spectrum of the closest Hll region, as SN 1995N sits in a irregular galaxy whose gradient is difficult to 
determine. 

^N2 metallicities obtained from 03N2 values converted to the N2 scale using the relation in|Kewley & Ellis^j2008). 


Note. — We assumed the metallicity ^adients measured by P04 for the hosts of eight SN IMs, and we checked them against our 
measurements, when available (see Fig.^. In the last column, “SN” indicates that the SN is still bright in the optical. All the NCR 
values are from H14. For SNe 2002kg andz003gm there are also metallicity measurements by |Maund et al.|(2006) , which are consistent 
with our results. 
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Table 6. Absolute (U/u, B, V, R/r/unf., I ji) peak magnitudes for the SNe Iln, Ibn and la-CSM of our sample. 


SN 

Type 

U j 'Umax 

^max 

4^ax 

i7/r/«n/.„ax 

/ ^max 



(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

1987F 

IIn-88Z 



< -17.38 

< -19.08 



1988Z 

Iln-88Z 



-18.40 

-18.68 

-18.97 


1994W 

lIn-P(94W) 



-17.20 

-17.33 

-17.31 


1994Y 

lIn-L(98S) 



-18.54 

-18.53 

-18.52 


1994ak 

lln 





< -16.82 


1995G 

Iln-88Z 

< 

-18.37 

< -17.97 

< -18.25 

< -18.50 

< -18.64 

1995N 

Iln-88Z 

< 

-13.68 

< -13.16 

< -13.80 

< -15.26 

< -13.99 

1997ab 

Iln-88Z 



< -18.82 




1998S 

lIn-L(98S) 



-18.31 

-18.37 

-18.47 

-18.46 

1999gb 

lln 





< -18.49 


2002A 

Iln 





< -16.22 


2002ea 

lln 





< -16.34 


2003ke 

lln 





< -16.44 


20031O 

lln 





< -16.66 


2003ma 

Iln-88Z 




-20.74 

-20.72 

-20.83 

2004F 

lln 





< -16.59 


2005bx 

lIn-L(98S) 



< -18.54 

< -18.06 

< -18.48 

< -18.67 

2005cp 

Iln-88Z 



-18.02 

-18.22 

-18.36 

-18.69 

2005db 

lln 



-15.72 

-16.48 

-17.08 

-17.65 

2005gl 

lIn-L(98S) 





-17.17 


2005ip 

IIn-88Z 

< 

-16.59 

< -16.87 

< -17.23 

< -17.38 

< -17.36 

2006aa 

lln 


-17.09 

-17.17 

-17.28 

-17.34 

-17.25 

2006am 

lln 





< -14.31 


2006bo 

IIn-P(94W) 

< 

-15.72 

< -16.29 

< -16.58 

< -16.70 

< -16.75 

2006db 

lln 





< -18.10 


2006jd 

Iln-88Z 

< 

-16.04 

< -16.42 

< -16.80 

< -17.07 

< -16.94 

2006qq 

Iln-88Z 

< 

-17.29 

< -17.62 

< -17.91 

-17.99 

-17.96 

2007n 

Iln-88Z 

< 

-17.67 

< -17.26 

< -17.66 

< -17.97 

< -18.02 

2007va 

lln 





< -18.47 

< -18.46 

2008B 

lln 





< -18.12 


2008fq 

IIn-L(98S) 

< 

-16.93 

-17.47 

-18.00 

-18.15 

-18.21 

2009ip 

IM-^IIn 


-19.18 

-18.19 

-18.22 

-18.31 

-18.39 

2010jl 

Iln-88Z 

< 

-19.64 

< -19.14 

-19.61 

< -19.894 

-20.24 

201 Iht 

IIn-P(94W) 


-17.71 

-16.70 

-16.67 

-16.77 

-16.83 

2013L 

Iln-88Z 


-19.11 

-19.02 

-19.03 

-19.01 

-18.88 

2008J 

la-CSM 

< 

-17.55 

-18.66 

-19.56 

-20.08 

-20.21 

1999cq 

Ibn 





-19.53 


2002ao 

Ibn 



< -15.25 

< -15.42 

< -15.68 

< -15.73 

20051a 

Ibn/IIn 





< -17.07 


2006jc 

Ibn 

< 

-18.65 

< -17.73 

< -17.56 

< -18.03 

< -17.90 

2010al 

Ibn 

< 

-19.41 

-18.55 

-18.59 

-18.743 

-18.80 

LSQ12btw 

Ibn 




< -18.96 

< -17.76 

< -17.64 
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Table 7. Wind velocity and mass-loss rate for a subsample of CSI SNe (Iln, Ibn, la-CSM). 


SN 

Type 

(km s“^) 

M 

(M© yr-i) 






Reference 


1987F 

IIn-88Z 

150 

0.01 

Wegner & SwansonI 

996 

IChugail) 

1991 

1 

1988Z 

IIn-88Z 

<200 

0.0007-0.015 

StathaKis & Sad 

ler 

m 


Lliugai & 

Uanziger |iyy4| 

1994W 

IIn-P(94W) 

1000 

0.3 

Sollerman et al. 

11 

59snch 

ugai et al. 

l2004bl "-' 

1995G 

IIn-88Z 

1000 

0.1 

Pastorello et al. 

[2002 

'.Itviewe et al.l 

20L 

1 dChugai & DanzigerH2003| 

1995N 

IIn-88Z 

<500 

0.0002 

l^ransson et al.l 

2002 

Izampieri et al 


a 

1997ab 

IIn-88Z 

90 

0.01 

Salamanca e 

t al 

IJiyyi 

0 

1998S 

IIn-L(98S) 

30-100 

0.0001 

lable y in K. 

lew 

e et al. |2012j, Anupama et al. |200lj ' 

2003ma 

IIn-88Z 

110 

0.13 

Rest et al. 1 

01 

r 






2005bx 

IIn-L(98S) 

813±133 

0.037±0.019 

Riewe et a . 

2(1 

12 






2005cp 

IIn-88Z 

eioitiio 

0.026±0.005 

Kiewe et al. 

Tun 






2005db 

Iln 

958±113 

0.057±0.024 

Kiewe et al. 

2012 






2005gl 

IIn-L(98S) 

420 

0.03 

Ual-Yam & 

Leona] 

LdJJ 

2009 

1 



2005ip 

IIn-88Z 

100-200 

0.0002-0.0004(0.0023) 

Smitti et al.l 

ZUUybUladdia et al.NzUlia 

Moriya et al.||2014| (lyr) 

2006aa 

Iln 

600 

0.005-0.016(0.8) 

laddia et al. 

2U13a .IMonva etal.||2Ul4 


2006bo 

IIn-P(94W) 

600 

<0.023-0.034 

laddia et al. 





2006jd 

IIn-88Z 

<160 

<0.007-0.032(0.021) 

S12,|'lhddia 


2013a|,| 

VIoriya et al.l|2014| 

2006qq 

IIn-88Z 

<200 

<0.003-0.007 

Taddia et al. 


a| 





2007rt 

IIn-88Z 

128 

0.056 

Irundle et a 

1)20091 





2008fq 

IIn-L(98S) 

500 

0.0011 

laddia et al. 

fzUljal 





2009ip 

Iln 

200 

0.009 

Margutti et 

il. 

7mA , 

Moriya etaljj201^ 

2010jl 

IIn-88Z 

100 

0.1 

Pransson et al. 

2014 


201 Iht 

IIn-P(94W) 

500-600 

0.01(0.03) 

Mauertian et al. 

2013 

D|,|Moriya et al.|j2014| 

2008J 

la-CSM 

50-100 

0.003 

laddia et ai. 

m 

■m— 

1999cq 

Ibn 

2000 


Matheson et al. 


00 





2002ao 

Ibn 

900-1500 


Pastorello et al. 

2DDK 

a 




20051a 

Ibn/IIn 

1500 


Pastorello et al. 

20011 

5 




2006jc 

Ibn 

600 


Anunama et al. 






2010al 

Ibn 

1000 


Pastorello et al. 

2015 

b 




LSQ12btw 

Ibn 

3200 


Pastorello et at. 

2015 

a 





jpassia et al.[|200l| report a mass-loss rate of 2 X 10 ® Mq yr ^ for the outer CSM and 3 X 10 ^ Mq yr ^ 
by[Anupam^^T.|^2001| is in the middle of this range. 


for the inner CSM. The value reported 
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★ 

ig54J/2002kg 

♦ 

1961V 

■ 

1987F 

■ 

1988Z 

• 

1994W 

♦ 

1995G 

■ 

1995N 

♦ 

1996L 

♦ 

1997bs 

• 

1998S 

★ 

1999bw 

■ 

1999cq 

♦ 

2000ch 

♦ 

2001ac 

■ 

2002ao 

★ 

2002bu 

★ 

2003gm 

★ 

OT2005 

◄ 

2005bx 

♦ 

2005cp 

• 

2005dp 

♦ 

2005gl 

* 

2005kj 

■ 

20051a 

■ 

2006aa 

• 

2006bo 

♦ 

2006fp 

• 

2006jc 

◄ 

2006qq 

★ 

2007SV 

• 

2008J 

♦ 

2008S 

◄ 

2008fq 

♦ 

OT2009 

• 

2010al 

♦ 

2010dn 

◄ 

PTFlliqb 

♦ 

2013L 

★ 

VI 



Figure 4 BPT diagram for each H II region that we spectroscopically observed. Above the dashed line ([Kauffmann et al.|2003 


Sect. 1^ the line fluxes are AGN contaminated and we rejected those spectra from our metallicity study. Most of the spectra fall 
in the star-forming part of the BPT diagram, below the dashed line. The symbols marked with black edges indicate that the ratio 
between [O III A5007] and H/3 is an upper limit. All the line ratios are reported in Table 
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Figure 5 Metallicity gradients of 10 SN Iln host galaxies observed with the NOT+ALFOSC. Symbols and lines as in the bottom 
panel of Fig.[^ Black open squares with arrows correspond to upper limits. The gradients of the hosts of SNe 2006jd and 2005ip are 
not shown as we could not measure the flux from bright H ll regions so we used literature data. For the host of SN 1995N, whose 
gradient is not shown, we only have one measurement and we did not assume any gradient, as it is an irregular and interacting 
galaxy. 
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Figure 6 Metallicity gradients of 5 SN Ibn and 1 SN la-CSM host galaxies observed with the NOT+ALFOSC. Symbols and lines 
as in the bottom panel of Fig. 
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Figure? Metallicity gradients of 10 SN IM host galaxies observed with the NOT+ALFOSC. Symbols and lines as in the bottom 
panel of Fig. We label the gradient “P04” when we use the one measured by [Pilyugin et al. ( 2004| l, and “<P04>” when we 
assumed the average P04 gradient. The data of the hosts of VI and OT2005 are not shown as for them we do not determined or 
adopt a gradient, but we only measured the metallicity of one H II region close the SN location. 
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Figure 11 Absolute peak magnitudes for SNe Iln and Ibn against the metallicity at the SN position. Triangles are upper limits. 
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Figure 12 (Top-left panel) Mass-loss rates of SNe Iln versus local metallicity measurements. Triangles are upper limits. Above 
10“^ Mq yr“^, the mass loss rates appear to be higher at higher metallicity. The power-law (PL) dependence on metallicity 
(a = 0.69 (e.g., [\^[20TT |)) of the line-driven mass-loss rates is shown with a black segment and is consistent with the data. 
SNe Iln-L (98S-like) typically show lower mass-loss rates than long-lasting SN Iln (88Z-like). (Top-right panel) Mass-loss rate 
CDFs for metal-poor and metal-rich SNe Iln, shown to better highlight that at higher metallicities the mass loss rates tend to be 
higher. Two different metallicity cuts are shown, at log(0/H)-Fl2 = 8.3 and 8.5. (Bottom-left panel) Wind velocities of SNe Iln 
and SNe Ibn versus local metallicities. Triangles are upper limits. The PL dependence on metallicity (a = 0.12 (e.g., Kudritzki 
2002|l) of the line-driven wind velocities rates is shown with a black segment, and is consistent with the data. SNe Iln-L (98S-like) 


typically show higher wind velocities than long-lasting SN Iln (88Z-like). (Bottom-rightpanel) Wind-velocity CDFs for metal-poor 
and metal-rich SNe Iln, shown to better highlight that at higher metallicities the wind velocities tend to be slightly higher. Two 
different metallicity cuts are shown, at log(0/H)+12 = 8.3 and 8.5. 
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Figur e 13 (Left-handpanel) Metallicity CDFs for CSI SNe showing and not showing MIR dust emission at late epochs (Fox et ah] 


2011 1 . These CDFs do not show any statistically significant difference. (Right-hand panel) CDFs of the peak absolute rjRfunf. 
band magnitudes for the same CSI SNe showing and not showing MIR dust emission at late epochs. These CDFs show that SNe 
with MIR dust emission tend to be brighter at peak. 
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Figure 14 (Toppanel) SN IM peak absolute magnitudes (from Smith et al.|201 la| l versus local metallicity. (Centralpanel) Absolute 
magnitude of SN IM progenitors (from Smith et al.|201 la| versus local metallicities. Triangles are lower limits. With the exception 
of SN 1961V and of those events where we have limits poorer than —12 mag, at lower metallicity the luminosity seems to be higher. 
(Bottom panel) The difference between each peak outburst and the corresponding progenitor magnitude versus metallicity is shown. 
Triangles are lower limits. The data suggest larger AM at larger metallicity. Red corresponds to R/unf. bands, blue to B band, 
green to V band, magenta to I band, brown to NIR, orange to R/V. 
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Table 3. Deprqjected and normalized distance from the host center and line ratios for the H II regions with observed spectrum. 


SN r/R25 

logioC F([N1I] A6584)/F(Ha)) 

logio(F([Om] A5007)/F(H/3)) 

2006aa 0.018 

-0.436 

-0.550* 

0.304 

-0.347 

-0.464* 

0.366 

-0.434 

-0.741* 

2002ao 0.067 

-0.913 

0.121 

0.073 

-0.818 

0.438 

0.204 

-1.012 

0.297 

0.315 

-1.016 

0.162 

20051a 0.022 

-0.494 

-0.178 

0.163 

-0.782 

0.017 

0.321 

-0.842 

0.161 

0.504 

-0.800 

0.119 

1988Z 0.076 

-0.714 

-0.034 

0.151 

-0.723 

0.021 

0.621 

-1.426* 

0.278* 

1999cq 0.005 

-0.283 

-0.564* 

0.351 

-0.400 

-1.185* 

0.389 

-0.524 

-0.363* 

0.524 

-0.449 

-0.083* 

0.573 

-0.535 

-2.721* 

0.654 

-0.570 

-0.319* 

1987F 0.016 

-0.449 

-0.814* 

0.201 

-0.677 

-1.662* 

0.224 

-0.456 

-0.544* 

0.348 

-0.517 

-0.181* 

0.422 

-0.522 

-0.292 

0.537 

-0.674 

-0.037 

1995N 

-0.938 

0.299 

2006jc 0.194 

-1.051 

0.363 

0.349 

-1.053 

0.322 

1994W 0.000** 

-0.246 

-0.371 

0.054 

-0.450 

-0.636 

0.090 

-0.558 

-0.916* 

0.142 

-0.639 

-1.002* 

0.173 

-0.502 

-1.028* 

0.209 

-0.509 

-0.834* 

0.328 

-0.546 

-3.062 

0.394 

-0.476 

-0.167 

1998S 0.397 

-0.534 

-0.808* 

0.401 

-0.571 

-0.485* 

0.407 

-0.539 

-0.477* 

0.438 

-0.588 

-0.662* 

0.458 

-0.515 

-0.394* 

0.517 

-0.515 

-0.394* 

0.531 

-0.330 

-0.400* 

0.556 

-0.470 

-0.570* 

0.598 

-0.367 

-0.520* 

0.679 

-0.367 

-0.520* 

0.725 

-0.555 

-0.213* 

0.756 

-0.497 

-0.292* 

2010al 0.064 

-0.498 

-0.470 

0.138 

-0.530 

-0.505 

0.150 

-0.560 

-0.830 

0.237 

-0.524 

-0.475 

0.392 

-0.656 

0.234 

0.536 

-0.660 

-0.044 

2005dp 0.000** 

0.265 

1.019 

0.283 

-0.724 

-0.883* 

0.395 

-0.577 

-1.101* 

0.532 

-0.515 

-1.181* 

0.906 

-0.397 

-2.315* 

2006bo 0.044 

-0.474 

0.023 

0.431 

-1.636* 

0.128 

0.654 

-1.478* 

0.497 

2006bo 0.029 

-0.397 

-0.228* 

0.114 

-0.594 

-0.573* 

0.190 

-0.548 

-1.216* 

0.214 

-0.336 

-0.417* 

0.249 

-0.498 

-0.594 

0.262 

-0.468 

-0.289* 

0.314 

-0.517 

-0.175 

0.352 

-0.581 

-0.067 

0.529 

-0.839 

-0.059 

2005gl 0.000** 

0.362 

0.367 

0.388 

-0.540 

-0.942* 

0.694 

-0.488 

-1.711* 

2005cp 0.031 

-0.535 

-0.898 
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Table 3 (cont’d) 


SN 

r/R25 

logiot F([N II] A6584)/F(Ha) ) 

logio(F([OiII] A5007)/F(H/3)) 


0.109 

-0.608 

-0.794* 


0.166 

-0.607 

-0.551 

1961V 

0.281 

-0.528 

-0.180 


0.428 

-0.467 

-0.245 


0.452 

-0.641 

0.102 

OT2009 

0.000 

-1.159 

0.192 


0.257 

-1.327 

0.549 


0.327 

-1.027 

0.378 


0.530 

-1.021 

0.239 


0.634 

-1.068 

0.361 

I995G 

0.000 

-0.231 

-0.323* 


0.404 

-0.720 

0.014 


0.476 

-0.678 

-0.158 


0.911 

-0.720 

0.170 


0.958 

-0.535 

-0.047 

2006fp 

0.480 

-1.191 

0.442 

2008S 

0.249 

-0.485 

-0.333* 

2003gm 

0.310 

-0.769 

-0.360 


0.340 

-0.612 

-0.381 

OT2005 

SN 

-2.033 

0.430 

1999bw 

0.349 

-0.661 

-0.116 


0.349 

-0.592 

-0.244 


0.350 

-0.659 

-0.054 


0.352 

-0.586 

0.149 


0.353 

-0.699 

0.124 


0.357 

-0.580 

-0.214 


0.394 

-0.650 

-0.046 


0.542 

-0.709 

0.123 

2002bu 

0.482 

-0.821 

0.487 


0.609 

-0.873 

-0.943* 

VI 

SN 

-1.937 

0.786 

1954J 

0.143 

-0.699 

0.123* 


0.222 

-0.964 

0.110 

2002kg 

0.042 

-0.708 

-0.514 


0.051 

-0.563 

-0.279 


0.065 

-0.382 

-0.190* 


0.070 

-0.893 

0.118 


0.102 

-0.948 

0.069 


0.131 

-0.787 

-0.324 


0.143 

-0.642 

-0.185 


0.174 

-0.832 

-2.526* 

2007SV 

SN 

-1.644 

0.685 

2000ch 

0.047 

-0.880 

0.189 


1.039 

-1.192 

0.020 


1.077 

-0.984 

0.291 

1997bs 

0.435 

-0.544 

-1.015 


0.447 

-0.528 

-0.902 


0.516** 

-0.544 

0.649 

200lac 

0.000** 

-0.248 

-0.250 


0.225** 

0.024 

0.010 


0.346 

-0.365 

-0.436 


0.968 

-0.353 

-0.215 

200lac 

0.000 

-0.504 

-0.792 


0.287 

-0.566 

-0.586* 


0.292 

-0.645 

-0.696* 


0.364 

-0.573 

-2.675* 


0.384 

-0.509 

-0.088* 


0.467 

-0.418 

-0.153 

2013L 

0.000 

-0.455 

-0.248 

I996L 

0.000** 

-0.354 

0.378 

2005kj 

0.000** 

-0.321 

-0.244 

2006qq 

0.000 

-0.369 

-0.201* 

2006qq 

0.000 

-0.346 

-0.631 

2005bx 

0.000 

-0.368 

-0.380 

PTFlliqb 

0.000** 

-0.281 

0.000 


* Upper limit 

***Excluded from the metallicity analysis after spectral classification with the BPT diagram. 



